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Introduction 


We  are  developing  hyperpolarized  MRI  allowing  for  ultrafast  (potentially  sub-second)  molecular  imaging  of 
breast  cancer  with  specificity  similar  to  that  of  FDG-PET,  but  significantly  lower  cost  than  conventional  MRI  due 
to  much  higher  patient  throughput/much  faster  exams  and  the  use  of  low  field  =  low  cost  MRI.  We  focus  our 
research  effort  on  the  high-risk  and  critical  challenges  that  must  be  solved  to  enable  clinical  implementation  of 
hyperpolarized  metabolic  imaging  of  breast  cancer. 

To  enable  efficient  imaging  of  these  agents,  this  proposal  explores  the  fundamental  property  of 
hyperpolarized  contrast  agents,  namely  large  nuclear  spin  polarization,  which  is  independent  of  the  applied 
magnetic  field  used  for  MRI.  This  key  feature  of  these  agents  allows  for  sub-second  imaging  with  smaller 
(lower  magnetic  field),  less  costly  scanners,  while  still  achieving  sensitivity  comparable  to,  or  greater  than  that 
of  conventional  high  field  MRI  scanners. 

Specifically,  the  research  efforts  during  Year  1  have  focused  on  two  specific  aims  as  described  in 
Statement  of  Work: 

1)  Preparation  of  pure  hyperpolarized  contrast  agents  approved  for  injection  in  humans 

l.a.  Y1Q1-Y1Q3.  Design  of  2nd  generation  129Xe  polarizer.  The  design  took  advantage  of  the  design  and 
development  process  of  the  1st  generation  XENA  device  ( 1 ,  2),  also  see  Appendix  3.  Xenon-129  polarizer  is 
needed  to  perform  Xenon-Induced  Polarization  of  13C-acetate  and  13C-pyruvate  that  will  serve  as 
hyperpolarized  contrast  agents  of  molecular  imaging  of  breast  cancer. 

I.b.  Y1Q4-Y2Q2.  Construction  of  the  open-source  2nd  generation  129Xe  polarizer  is  an  on-going  process,  which 
is  reported  here. 

2)  Development  of  low  field  (low  cost)  sub-second  in  vivo  imaging  suitable  for  preclinical  models  of 
breast  cancer  and  clinical  trials 

4.a.  Y1Q1-Y1Q2.  Low  field  MRI:  pre-clinical  development  was  done  based  on  our  existing  hardware  that  was 
developed  during  our  preliminary  work.  The  work  focused  on  developing  and  construction  of  multiple  RF  coils 
intended  for  optimal  direct  and  indirect  detection  of  hyperpolarized  contrast  agents  in  vivo. 

4.b.  Y1Q3-Y1Q4.  Low  field  MRI:  pre-clinical  optimization  was  conducted  to  optimize  developed  hardware  and 
MRI  sequences  for  specific  metabolic  contrast  agents  and  their  metabolites  in  vivo.  Specifically,  we  optimized 
the  detection  sensitivity  and  developed  a  theoretical  framework  for  future  clinical  scale  RF  coil  designs  as  well 
as  limits  of  detection  sensitivity  for  low-field  MRI. 
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Body 


Task  l.a.  129Xe  Polarizer:  design 

The  second-generation  xenon-129  polarizer  design  was  largely  based  on  the  first-generation  device  with 
significant  improvements  that  enhance  stability,  polarization  yield  of  the  contrast  agent,  automation,  ease  of 
operation/maintenance,  and  future  integration  with  xenon-induced  polarization  module  and  smaller  footprint. 
The  device  consists  of  the  following  main  components: 

1.  Gas-handling  manifold; 

2.  VHG-narrowed  200  W  laser  diode  array  (LDA)  module; 

3.  Variable  temperature  oven; 

4.  Electromagnet; 

5.  Low-field  NMR  spectrometer  with  in  situ  detection; 

6.  In  situ  infrared  (IR)  spectroscopy; 

7.  Micro-controller; 

8.  Integrating  frame  /  chassis; 

The  design  was  largely  performed  using  SolidWorks  CAD  software  available  to  our  team  free  of  charge  via  the 
core  facility  at  the  Imaging  Institute.  Provided  below  is  the  description  of  the  components’  design. 

1.  Gas-handling  manifold 

The  gas-handling  manifold  achieves  oxygen-free  operation  by  always  maintaining  high-pressure  (i.e.  above 
atmospheric)  conditions  upstream  of  the  cell,  as  shown  in  Fig.  1 .  The  manifold  is  mounted  on  the  aluminum 
plate  that  is  precision  cut  as  a  part  of  the  component  8,  the  integrating  frame. 

XEUS-MRKII  POLARIZER  MANIFOLD  SCHEMATIC 
(OPTIONAL:  NITROGEN  For  Purge  Cycles) 
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Figure  1.  Gas  Loading  Manifold  Design.  By  using  a  pre-mixed  Xe/N2  gas  cylinder,  we  have  significantly 
reduced  the  number  of  pneumatic  valves  used  from  past  designs — while  also  mitigating  the  purge  cycle 
loading  procedure,  which  vastly  reduces  the  time  required  to  load  the  Optical  Pumping  (OP)-Cell  in  preparation 
for  the  next  run  (increasing  the  overall  HP129Xe-generating  efficiency  of  the  device).  This  approach  also 
reduces  the  risk  of  02  contamination  by  always  leaving  the  loading  lines  pressurized  to  the  desired  loading 
pressure,  which  are  programmed  via  the  microcontroller  to  re-pressurize  if  a  drop  in  pressure  is  detected.  The 
OP  cell  at  the  center  of  the  polarizer  is  a  500  cc  Pyrex  cylinder  with  two  stopcock-sealed  valves  to  be 
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controlled  by  rotary  actuators.  Many  such  cells  have  now  been  delivered  and  may  now  be  purchased  as  a 
stock  part  from  our  supplier  (Midrivers  Glass)  guaranteed  to  fit  our  oven  design  (despite  the  otherwise 
inherently-custom  nature  of  such  glass  pieces).  Prior  to  gas  loading  and  install,  the  cells  are  cleaned,  coated 
with  a  siliconizing  agent  (Surfasil),  loaded  with  -200  mg  Rb  in  an  inert-atmosphere  glove  box,  and  the  Rb  is 
distributed  as  a  thin  film  throughout  the  cell. 

2.  VHG-narrowed  200  W  LDA  module 

The  200  W  laser  diode  array  (LDA)  module  consists  of  an  LDA  (internally  narrowed  via  volume  holographic 
grating  (VHG)  technology)  purchased  from  a  commercial  vendor  (QPC),  which  was  approved  in  the  original 
budget.  We  worked  closely  with  the  company  to  ensure  that  the  laser  module  has  an  integrated  optical  train 
that  delivers  polarized  light  and  expands  it  into  a  2”  OD  beam  that  homogenously  irradiates  the  2”  optical 
pumping  cell,  Figs.  1-2.  The  module  also  contains  two  heat  sinks  providing  efficient  heat  dissipation.  The  heat 
sink  fans  are  gated  by  the  micro-controller  to  suppress  the  background  signal  for  high  quality  QA  by  in  situ 
NMR  spectroscopy.  The  laser  temperature,  which  must  be  maintained  constant  to  ensure  operation  at  a 
desired  output  wavelength  for  optimal  absorption  by  the  Rb  vapor,  is  controlled  using  a  1,400  W  water  chiller. 
Because  laser  performance  is  a  key  to  overall  polarization  success  of  xenon-129  polarizer  and  because  it  is 
the  first-of-its-kind  integrated  device  made  by  this  vendor,  the  SIUC  team  (Prof.  Boyd  Goodson  lab)  tested  the 
performance  of  this  laser  and  the  results  of  this  test  are  summarized  in  Appendix  2. 

200W  VHG-Narrowed  3D  Printed  Oven 

LDA  LASER  Optical  Pumping  Oven 


Figure  2.  200  W  VHG-Narrowed  LDA  laser  schematic  aligned  concentrically  with  the  optical  pumping  oven. 


3.  Variable-temperature  (VT)  oven 

The  VT  oven,  Fig.  3,  controls  the  surface  temperature  of  the  OP  cell  to  maintain  the  Rb  vapor  density  during 
spin-exchange  optical  pumping,  the  process  by  which  spin  angular  momentum  from  circularly-polarized 
photons  is  converted  into  Rb  electronic  spin  polarization,  and  ultimately,  to  129Xe  nuclear  spin  polarization.  The 
oven  was  designed  at  VUIIS  using  SolidWorks  CAD  software  and  3D  printed  in  polycarbonate  material  in  the 
laboratory  of  our  collaborator  Prof.  Matthew  Rosen  (Harvard  University).  The  oven  also  features  a 
thermoelectric  cooler  (TEC)/heater  that  maintains  oven  temperature  by  controlling  the  temperature  of 
flowing/re-circulating  air  pushed  by  an  integrated  blower.  The  resulting  design  does  not  require  any 
chilling/cooling  source  besides  electric  power.  That  provides  a  clear  advantage  over  the  previous  polarizer 
design  ( 1 ,  3).  The  TEC  module  is  controlled  by  a  micro-controller  for  future  complete  automation  during  Year  2. 
The  front  on  the  oven  has  an  optical  flat  allowing  efficient  transmission  of  the  laser  beam,  while  the  back  of  the 
oven  has  a  reflecting  mirror  to  maximize  the  use  of  the  laser  flux  (by  providing  a  second  pass  of  the  beam 
through  the  OP  cell)  and  improve  the  homogeneity  of  laser  illumination  along  the  OP  the  cell.  The  Oven  also 
houses  an  RF  coil  for  in  situ  NMR  and  a  fiber  optical  cable  in  the  back  of  the  oven  (not  shown)  for  in  situ  IR 
spectroscopy.  The  NMR  reports  on  the  quality  of  the  xenon-129  polarization,  while  the  IR  spectroscopy  reports 
on  the  Rb  polarization  and  the  laser  beam  absorption. 
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Figure  3.  Variable  temperature  (VT)  oven  housing  the  0.5L  (500  cc)  optical  pumping  (OP)  cell, 


3D  (Poly.  Carb.)  PRINTED  OPTICAL  PUMPING  OVEN  W/ 
THERMO-ELECTRIC  COOLER 


3D  CAD  SOFTWARE  DESIEGNED 
AT  VANDERBILT  -  NASHVILLE,  TN 


52  Hours  Later  ^ 


POLYCARBONATE  3D  PRINTED 
AT  HARVARD  -  BOSTON.  MA 


2”  OD.  The 


design  and  production  using  3D  printer  is  shown  on  the  left  and  the  exploded  view  is  shown  on  the  right. 


4.  Electromagnet 

The  electromagnet  utilizes  a  more  efficient  Barker  coil  design,  Fig.  4.  The  design  provides  a  homogenous 
magnetic  field  B0  for  the  optical  polarization  process  and  in  situ  NMR  detection. 


Figure.  4.  Electromagnet  design  using  four  individual  23.6”  OD  coils  comprising  two  217  (outer)  and  two  96 
(inner)  turn  coils.  This  Barker  electromagnet  coil  design  provides  a  sufficiently  homogenous  magnetic  field  that 
envelops  the  region  of  the  OP  cell  while  also  embodying  a  more  compact  design  (smaller  footprint)  compared 
to  other  coil  configurations.  The  individual  coils  made  by  Acutran  LP  (Fombell,  PA)  and  a  custom-designed 
mechanism  allows  for  primitive  shimming  of  the  magnetic  field  by  moving  the  individual  magnet  coils  with 
respect  to  each  in  four  corners  of  the  frame.  The  magnet  coils  are  designed  around  the  plane  that  houses  VT 
oven. 


5.  Low-field  NMR  spectrometer  with  in  situ  detection 

Kea2  dual-channel  NMR  spectrometer,  Magritek,  New  Zealand,  with  built-in  1  W  (per  channel)  RF  amplifier  is 
used  for  in  situ  NMR  detection  capability.  The  spectrometer  utilized  an  RF  coil  operating  at  84  kHz  and  very 
recently  was  re-tuned  to  47  kHz  (work  in  progress).  The  RF  coil  is  housed  next  to  the  OP  cell  and  reports 
onthe  xenon-129  polarization.  The  advantage  of  lower-frequency  detection  is  lower  operating  magnetic  field 
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(which  requires  less  energy  and  generates  less  waste  heat).  The  NMR  spectrometer  can  be  triggered  by  the 
micro-controller. 

6.  In  situ  infrared  (IR)  spectroscopy 

In  situ  IR  spectroscopy  is  enabled  with  a  HR4000  High-resolution  User-configured  Spectrometer,  Ocean 
Optics,  Dunedin,  FL.  This  high-resolution  IR  spectrometer  allows  for  monitoring  the  laser  output  spectral  profile 
(including  centroid  wavelength,  linewidth,  and  other  features),  as  well  as  the  Rb  absorption  spectral  profile 
(including  the  fraction  of  laser  light  absorbed  by  the  OP  cell). 

7.  Micro-controller 

The  polarizer  contains  a  micro-controller  that  utilizes  the  same  design  as  previously  described  ( 1 ,  3),  except 
that  its  shell  was  designed  in  SolidWorks  CAD  software  for  better  integration  with  the  Frame  /  Chassis  (see 
below). 

8.  Integrating  Frame  /  Chassis 

The  frame  was  designed  using  SolidWorks  CAD  software  and  precision  cut  to  our  design  specifications  by  a 
private  vendor  (Minitec  Systems  LLC,  Victor,  NY).  The  frame  integrates  all  components  described  above  as  well  as 
other  components:  a  silent  air-compressor  for  actuating  pneumatic  valves,  a  vacuum  pump  (Edwards,  vacuum  level  <  10 
mTorr)  and  other  minor  electrical  and  ventilation  components,  Fig.  5.  The  frame  also  includes  the  mounting  plates  for  the 
gas-handling  manifold. 


Figure  5.  Polarizer  3D  schematic  rendering  made  during  the  design  stages  of  the  XeUS  polarizer. 
The  overall  design  of  xenon-129  polarizer  was  completed  during  Year  1  of  the  project. 

Task  l.b.  129Xe  Polarizer:  construction 


1.  Completed  construction  Work 

While  the  polarizer  will  be  completed  during  the  second  year  of  the  project,  we  have  made  significant  progress 
with  its  construction.  Nearly  all  components  have  been  installed,  Fig.  6,  allowing  for  safe  device  operation  in 
the  manual  polarization  mode.  The  manual  mode  of  operation  does  not  currently  allow  recharging  the  cell  with 
new  batches  of  gas,  but  allows  for  device  optimization,  i.e.  mapping  optimal  conditions  for  achieving  maximum 
xenon-129  polarization.  Two  xenon  density  conditions  were  tested  during  Year  1:  xenon  :  nitrogen  12.5%  : 
87.5%  and  50%  :  50%  with  total  in  cell  pressure  of  2000  Torr.  These  cells  were  prepared  in  the  lab  of  co-PI 
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Prof.  Boyd  Goodson  (SIU  site).  The  preliminary  optical  pumping  results,  Fig.  7,  are  very  encouraging  indeed 
with  record  xenon-129  polarization  levels,  Fig.  8.  The  levels  of  xenon-129  polarization  under  these  two 
conditions,  Fig.  8c,  are  approaching  the  theoretical  maximum  of  100%.  The  in  situ  IR  measurements  of  Rb 
polarization,  Fig.  9,  are  in  accord  with  the  high  xenon-129  polarization  levels. 

It  should  be  noted  that  nearly  all  components  have  been  assembled,  although  integration  was  not  yet 
completed.  The  device  integration  and  automation  will  be  completed  after  optimal  operating  conditions  are  be 
determined.  Moreover,  we  have  already  demonstrated  record  levels  of  xenon-129  polarization  in  the  partially 
completed  device. 


Figure  6.  Left:  design.  Right:  device  under  construction.  XeUS  Polarizer  during  testing  phase:  (a)  200  W  VHG 
frequency-narrowed  QPC  laser  with  optical  train  aligned  with  OP  cell  producing  2”  diameter  beam  @794.76 
nm  with  a  bandwidth  of  ~0.3  nm.  (b)  3D  Printed  Optical  Pumping  Oven,  (c)  Gas  loading  manifold  components, 
(d)  Microcontroller  box,  responsible  for  controlling  all  pneumatic  valves  via  GUI  interface  allowing  manual  and 
automated  operation,  (e)  PSUs  for  Laser  and  Magnetic  Coils,  (f)  Kea2  low-field  NMR  Spectrometer,  (g)  1400 
Watt  Water  Chiller  for  laser  cooling,  (h)  Vacuum  System  and  Air  Compressor,  (i)  B0  Coil  with  Barker 
configuration. 
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Figure  7.  TEC  oven  during  optical  pumping  with  laser  light  illuminating  the  cell.  The  RF  coil  below  the  OP  cell 
samples  xenon-129  polarization  in  situ,  Fig.  8. 


Thermal 


Hyperpolarized  129Xe 


LDA  Output  =  1 50  Watts 
OP  Cells  Backfilled  with 


83  83.5  84  84.5  8  5  8  3  83.5  84  84.5  8  5  0  200  400  600  80010001200 

■•H  Larmor  Frequency  (kHz)  129Xe  Larmor  Frequency  (kHz)  Xe  In  Cell  Density  (Torr) 

:igure  8.  /n  situ  NMR  spectroscopy  at  84  kHz.  (a)  ^  reference  (Thermal)  spectrum  obtained  from  water, 

pH]  =  1 10.7  M,  doped  with  10  mM  CuS04  inside  0.50  L  OP  cell.  Spectrum  was  acquired  at  84  kHz  Larmor 
frequency,  with  100,000  averages  using  ^  thermal  polarization  P  =  6.8x1  O'9,  (b)  129Xe  spectrum  (nat.  ab.  = 
26.44%)  of  1000  torr  Xe,  [129Xe]  =  14.3  mM,  backfilled  with  N2  to  a  total  cell  pressure  of  2000  torr  acquired  at 
84  kHz  129Xe  Larmor  frequency  with  1  scan.  Both  spectra  were  taken  using  the  same  circuit-detection  coil  with 
the  same  B1  tipping  angle,  (c)  Bar  graph  representing  %Pxe  at  250  and  1 000  torr  with  -85%  and  74%  xenon 
polarization  respectively. 
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PXe  build-up  during  SEOP  at  62  QC;  the  red  curve  is  subsequent  129Xe  decay  after  the  laser  was  turned  off  and 
the  cell  was  cooled.  Both  datasets  were  fit  to  simple  exponentials  (curves);  red  data  points  were  obtained  with 
rf  pulses  possessing  half  the  tipping  angle.  Red  data  points  obtained  before  the  cell  reached  room  temperature 
(grey  highlight)  were  not  included  in  the  fit.  Losses  from  rf  pulses  were  ignored  during  the  fits,  (b)  In  situ 
monitoring  of  global  Rb  electron  spin  polarization.  Examples  of  near-IR  spectra  of  the  pump  laser  transmitted 
through  the  cell  acquired  under  the  following  conditions:  1000  Torr  Xe,  1000  Torr  N2;  room  temperature  (“cold 
cell”);  62  °C  and  B0=7.13  mT  at  the  beginning  of  SEOP  (“t=0”);  62  °C  and  B0=7.13  mT  at  steady-state  for  PXe 
(“t=SS”,  see  adjacent  figure);  and  62  °C  and  B0= 0  mT  at  steady-state.  The  large  changes  in  transmitted  laser 
intensity  when  the  field  is  cycled  indicate  high  Rb  electronic  spin  polarization  (see  inset). 

2.  Work  in  Progress  on  xenon-129  polarizer 

The  polarizer  construction  work  is  still  ongoing  from  the  perspective  of  identifying  the  optimal  operating 
conditions,  completing  device  automation  using  the  constructed  micro-controller  box,  and  testing  the 
constructed  gas-handling  manifold.  We  will  utilize  the  same  Rb  filter/getter  as  in  the  first-generation  device.  In 
case  if  it  will  be  insufficient,  a  cold  trap  (t<0  °C)  to  catch  residual  Rb  will  be  implemented  using  a  TEC-cooled 
trap. 


it 


Task  4.a.  Low  field  MRI:  pre-clinical  development 


We  have  developed  and  constructed  low-field  preclinical  imaging  system  for  mouse  imaging  in 
collaboration  with  Magritek,  New  Zealand.  While  the  imager  was  developed  prior  the  beginning  of  this  project, 
the  development  of  dedicated  RF  coil  for  low-field  MRI  was  essential  for  the  progress  of  preclinical  optimization 
from  the  perspective  of  optimization  of  signal-to-noise  ratio,  which  will  ultimately  determined  the  limits  of  spatial 
and  temporal  resolution.  Fig.  10  below  provides  the  summary  of  our  pre-clinical  MRI  system. 

A  L  R 


vrm 


Figure  10.  The  0.0475  T  MR  scanner  system.  A)  The  balanced  tuning/matching  circuit;  B)  Representation  of 
the  inner  basket  weave  detection  coil,  outer  Helmholtz  saddle  coil  for  ^CVH  excitation,  and  RF  shield 

support  tube;  C)  and  D)  views  of  the  0.0475  T  probe  design;  E)  realization  of  the  H-X  and  X-H  probes  with 
attached  tuning/matching  boxes;  F)  View  of  overall  MR  scanner  system  in  the  vertical  spectrometer 
configuration. 


Task  4.b.  Low  field  MRI:  pre-clinical  optimization 

The  optimization  was  largely  focused  around  optimization  of  detection  RF  coils.  MRI  acquisition  parameters 
were  also  determined.  Very  high  sensitivity  RF  coils  were  developed  as  well  as  the  theory  proving  that  low  field 
detection  of  hyperpolarized  contrast  agents  can  be  indeed  more  sensitive  than  that  at  high  magnetic  field  (4). 
Fig.  11  summarizes  the  imaging  work  at  low  and  high  field.  While  113%  of  high  field  SNR  was  theoretically 
possible,  the  limitations  of  our  imaging  system  limited  the  actual  RF  coil  performance  to  40%.  These  are 
important  lesson  that  will  be  taken  into  account  when  developing  human  scale  MRI  system. 
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B0  =  4.7  T  B0  =  0.0475  T 


X  (mm)  X  (mm) 

0  16  32  48  64  80  0  16  32  48  64  80 


Figure  11. 13C  and  1H  MRI  of  sodium  1-13C-acetate.  4.7  T  acquisition  used  Boltzmann  13C  P  =  4.06*1  O'6  and 
'H  P  =  1 .61*1  O'5,  and  0.0475  T  used  approximately  the  same  polarization  levels,  13C  P  =  4.70*1  O'6  and  'H  P  = 
1.05*1  O'5.  All  measurements  used  a  spherical  phantom:  for  ^  1.0  g  sodium  1 -13C-acetate  in  99.8%  D20  with 
2.8  mL  total  volume  and  for  13C  5.18  g  sodium  1-13C-acetate  in  99.8%  D20  with  17.5  ml.  total  volume.  All 
acquisition  and  processing  parameters  were  identical  except  13C  excitation  pulse  angle  a.  No  image 
extrapolation  or  zero  filling  was  used. 

We  note  the  sub-millimeter  resolution  level  achieved  in  proton  scans  of  sodium  1-13C-acetate.  For  future 
imaging  efforts  in  Year  2-3  in  preclinical  models  of  breast  cancer,  the  following  initial  imaging  parameters  will 
be  utilized. 

Imaging  sequence:  2D  GRE  without  slice  selection  using  direct  13C  detection 

TE/TR:  6.4  ms  / 14  ms 

Imaging  matrix:  32  x  32 

FOV:  80x80  mm2 

Spatial  Resolution:  2. 5x2. 5  mm2 

Total  Scan  Speed:  <0.5  s 

Excitation  pulse  angle:  15  degrees  for  a  single  image  acquisition  to  maximize  SNR  and  5  degree  excitation  RF 
pulse  for  dynamic  acquisitions  with  ~2  s  temporal  resolution. 

We  will  optimize  the  pulse  sequence  further  once  we  have  more  preclinical  experience  and  will  extend  the 
detection  to  indirect  proton  detection,  full  3D  acquisition  and  compressed  sensing  with  the  goals  of  optimizing 
concentration  limits  and  spatial/temporal  resolution. 

We  will  use  similar  or  identical  acquisition  parameters  on  pre-clinical  4.7  T  system  for  quantitative  comparison 
of  imaging  at  low  and  high  magnetic  fields. 
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Key  Research  Accomplishments 

•  Designed  second-generation  clinical  scale  xenon-129  polarizer 

•  Achieved  nearly  100%  theoretical  maximum  xenon-129  polarization  in  partially  constructed  polarizer 
demonstrating  its  future  success  and  utility  for  future  use  as  a  part  of  xenon-induced  polarization  setup  for 
polarization  of  13C  metabolic  contrast  agents 

•  The  in  situ  NMR  quality  assurance  resonance  frequency  of  developed  xenon-1 29  polarizer  was  decreased 
to  record  low  47  kHz  to  minimize  power  requirements  of  xenon-129  polarizer  and  minimize  its  thermal 
management 

•  Developed  novel  dual  channel  low-field  NMR  RF  coils  for  high  detection  sensitivity  of  hyperpolarized 
contrast  agents  in  pre-clinical  molecular  imaging  of  breast  cancer 

•  Achieved  40%  detection  sensitivity  of  MRI  at  0.0475  T  compared  to  that  of  4.7  T  pre-clinical  MRI  scanner 
for  direct  13C  detection  and  indirect  proton  detection.  More  that  100%  was  theoretically  possible,  but  not 
achieved  due  to  hardware  limitations 

•  Developed  quantitative  theory  of  MRI  detection  sensitivity  of  hyperpolarized  contrast  agents  in  low  field 
from  the  perspective  of  RF  coil  design.  This  will  pave  the  road  to  theoretical  RF  coil  design  for  future 
clinical  scale  system 

•  Developed  imaging  protocols  for  2D  in  vivo  molecular  preclinical  imaging  of  breast  cancer  in  mouse  model 
using  0.0475  T  MRI  scanner  with  optimized  RF  coils 

Reportable  Outcomes 

Oral  Presentations 

1 .  Chekmenev,  E.  Y.,  13  C  and 15  N  PH  IP  and  XI P  polarized  contrast  agents  for  molecular  imaging  of  cancer, 
COST  Action  TD1103  EuroHyperPol,  Jun.  29,  Dublin,  Ireland,  2012. 

2.  Chekmenev,  E.  Y.,  PHIP  polarizing  and  imaging  strategies  for  13C-lactate  and  15N-choline,  The  Third 
International  Workshop,  on  Metabolic  Imaging,  Jul.  27,  Philadelphia,  PA,  USA,  2012. 

3.  Chekmenev,  E.  Y.,  What  about  Choline,  IN  VIVO  MAGNETIC  RESONANCE  GORDON  RESEARCH 
CONFERENCE,  Aug.  3,  Colby  College,  Waterville,  ME,  USA,  2012. 

4.  Chekmenev,  E.  Y.  Molecular  Imaging  of  Breast  Cancer  using  Hyperpolarized  MRI  Contrast  Agents, 
Vanderbilt-lngram  Cancer  Center,  Breast  Cancer  Program  &  SPORE  Seminar,  Oct.  1,  Vanderbilt 
University  Medical  Center,  Nashville,  TN,  USA,  2012. 

5.  Chekmenev,  E.  Y.  Molecular  Imaging  of  Breast  Cancer  using  Hyperpolarized  MRI  Contrast  Agents, 
Vanderbilt  University  Institute  if  Imaging  Science  Seminar,  Jan.  1 1 ,  Vanderbilt  University  Medical  Center, 
Nashville,  TN,  USA,  2013. 

Conference  Abstracts 

1.  Chekmenev,  E.  Y.  (Jul.  27,  2012).  PHIP  polarizing  and  imaging  strategies  for  13C-lactate  and  15N-choline. 
The  Third  International  Workshop  on  Metabolic  Imaging,  Philadelphia,  PA,  USA. 

2.  Coffey,  A.  M.,  M.  Truong  and  E.  Y.  Chekmenev  (April  14-19,  2013).  Improving  imaging  sensitivity  of 
hyperpolarized  MR  at  low  magnetic  detection  field  strengths.  Experimental  NMR  Conference,  Asilomar, 
California. 

3.  Nikolaou,  P.  and  E.  Y.  Chekmenev  (2012).  An  ‘Open-Source’  129Xe  Polarizer  for  Clinical  Imaging  and 
Secondary  Hyperpolarization  of  13C  Tracers  for  Molecular  Imaging  of  Cancer  Vanderbilt  University 
Imaging  Institute  Annual  Retreat. 

4.  Nikolaou,  P.,  A.  Coffey,  L.  L.  Walkup,  B.  Gust,  H.  Newton,  I.  Muradyan,  M.  Rosen,  S.  Patz,  M.  J.  Barlow,  B. 
M.  Goodson  and  E.  Y.  Chekmenev  (April  14-19,  2013).  "XeUS"; — A  Second-Generation  ‘Open  Source’ 
Clinical-Scale  129Xe  Polarizer  for  Xenon-Induced  Polarization.  Experimental  NMR  Conference,  Asilomar, 
California. 

5.  Nikolaou,  P.,  A.  Coffey,  L.  L.  Walkup,  B.  Gust,  N.  Whiting,  H.  Newton,  S.  Barcus,  I.  Muradyan,  M. 
Dabaghyan,  G.  D.  Moroz,  M.  Rosen,  S.  Patz,  M.  J.  Barlow,  E.  Y.  Chekmenev  and  B.  M.  Goodson  (April  14- 
19,  2013).  High  (-30-90%)  129Xe  Hyperpolarization  at  High  Xe  Densities  Using  an  ‘Open-Source’ 

Polarizer  for  Clinical  and  Materials  MRS/MRI.  Experimental  NMR  Conference,  Asilomar,  California. 
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Manuscripts 

1.  R.  V.  Shchepin,  W.  Pham,  E.  Y.  Chekmenev,  Towards  Hyperpolarized  Phospholactate  Imaging: 
Dephosphorylation  and  Biodistribution  of  1-13C-phospholactate  In  Vivo.  Contrast  Media  Mol.  Imaging, 
submitted  (2013). 

2.  R.  G.  Abramson,  L.  R.  Arlinghaus,  J.  Weis,  i.  X.  L,  A.  N.  Dula,  E.  Y.  Chekmenev,  S.  Smith,  M.  I.  Miga,  V. 
Abramson,  T.  E.  Yankeelov,  Current  and  emerging  quantitative  magnetic  resonance  imaging  methods  for 
assessing  and  predicting  the  response  of  breast  cancer  to  neoadjuvant  therapy.  Breast  Cancer  Targets 
and  Therapy  4,  139  (2012). 

3.  P.  Nikolaou,  A.  Coffey,  L.  L.  Walkup,  B.  Gust,  N.  Whiting,  H.  Newton,  S.  Barcus,  I.  Muradyan,  M. 
Dabaghyan,  G.  D.  Moroz,  M.  Rosen,  S.  Patz,  M.  J.  Barlow,  E.  Y.  Chekmenev,  B.  M.  Goodson,  Near-unity 
nuclear  polarization  with  an  ‘open-source’  129Xe  hyperpolarizer  for  NMR  and  MRI.  Proc.  Natl.  Acad.  Sci. 

U.  S.  A.,  conditionally  accepted  (2013). 

Funding  applied  for  based  on  work  supported  by  this  award 

1 .  3D  Printer  for  Biomedical  Research  and  Engineering,  submitted  to  DoD  FY201 3  DEFENSE  UNIVERSITY 
RESEARCH  INSTRUMENTATION  PROGRAM  (DURIP),  Department  of  Army,  PI:  Chekmenev,  Requested 
amount:  $297,586.00,  Date  submitted:  2012-09-27. 

2.  MRI  Consortium:  Portable  Large  Scale  Low  Field  Xenon- 129  Polarizer/Sub-Second  MRI  Scanner, 
submitted  to  NSF  MRI  program,  Pis:  Chekmenev  (Vanderbilt),  Goodson  (Southern  Illinois  University), 
Requested  amount:  $  840,200,  Date  submitted:  2013-02. 

Conclusion 

Second-generation  xenon-129  polarizer  was  designed,  and  it  now  undergoes  the  construction  according  to 
the  DoD  approved  schedule.  The  second-generation  polarizer  includes  several  key  engineering  innovations. 
When  combined,  they  resulted  in  more  efficient  polarizer  operation.  For  example,  the  polarizer  is  a  stand-alone 
device  that  does  not  require  any  every  day  consumables  besides  electric  power.  So  far  this  has  not  been 
achieved  with  any  clinical  scale  xenon-129  polarizer,  so  it  definitely  represents  an  engineering  breakthrough, 
which  can  be  utilized  in  biomedicine  in  the  future.  Relevant  to  our  consortium  on  xenon-129  polarizers’ 
development  this  paves  the  road  to  the  next  generation  devices  that  will  be  completely  automated  and  can  be 
potentially  operated  as  push-button  devices  from  the  perspective  of  contrast  agent  production  of 
hyperpolarized  xenon-129.  Scientifically,  we  have  achieved  very  high  level  of  xenon-129  polarization,  which  we 
currently  attribute  to  more  careful  design  and  cell  temperature  management.  Future  work  on  optimization  will 
provide  better  understanding  during  QA  process  in  Year  2  of  the  project.  High  xenon-129  polarization  is  critical 
for  xenon-induced  polarization,  because  xenon-129  will  serve  as  a  source  of  nuclear  spin  polarization  and 
approaching  theoretical  maximum  with  xenon-129  would  potentially  allows  approaching  theoretical  maximum 
with  13C  hyperpolarized  contrast  agents. 

We  have  demonstrated  that  low-field  imaging  of  hyperpolarized  contrast  agent  can  be  as  sensitive  as  the 
high-field  detection  and  potentially  more  sensitive  (but  would  require  cryogenic  RF  coils  -  these  are  outside  the 
scope  of  the  project  and  our  current  expertise)  than  high  field  detection.  We  have  also  demonstrated  this 
concept  with  40%  detection  sensitivity,  which  is  in  a  good  agreement  due  to  known  hardware  limitations: 
magnet  bore,  non-ideal  wire  selection,  etc. 

The  strategy  of  testing  and  optimizing  RF  coil  sensitivity  of  low  field  at  preclinical  stage  paid  off,  because 
we  were  able  to  developed  a  quantitative  theory  for  SNR  of  hyperpolarized  contrast  agent  in  low  field.  The 
lessons  learned  about  geometry  and  components  interactions  will  be  useful  during  clinical  research  system 
design  and  construction  and  will  undoubtedly  allow  us  to  achieve  optimal  performance.  This  is  essential  as  the 
development  of  clinical  system  in  Years  3-5  is  one  of  the  ultimate  goals  of  the  project. 

No  changes  for  future  work  are  requested. 
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Appendices 

Appendix  1:  Statement  of  work,  Year  1. 


Milestones 

Year  1 
quarters 

1 

2 

3 

4 

l.a.  Open  Source  129Xe  Polarizer:  design 

X 

X 

X 

I.b.  Open  Source  129Xe  Polarizer:  construction 

X 

4. a.  Low  field  MRI:  pre-clinical  development 

X 

X 

4.b.  Low  field  MRI:  pre-clinical  optimization 

X 

X 

This  statement  of  work  is  divided  in  four  different  sections  matching  four  Specific  Aims  of  the  proposal. 
Relevant  to  the  Year  1  are  the  following: 

3)  Preparation  of  pure  hyperpolarized  contrast  agents  approved  for  injection  in  humans 

l.a.  Y1Q1-Y1Q3.  Design  of  open-source  2nd  generation  129Xe  polarizer.  The  design  will  take  advantage  of  the 
design  and  development  process  of  the  1st  generation  XENA  device.  More  details  are  provided  in  the 
description  of  the  sub-contract  to  Southern  Illinois  University. 

I.b.  Y1Q4-Y2Q2.  Construction  of  the  open-source  2nd  generation  129Xe  polarizer.  The  construction  will  require 
low  field  spectrometer  and  high  power  narrowed  tunable  laser  requested  in  the  original  budget. 

4)  Development  of  low  field  (low  cost)  sub-second  in  vivo  imaging  suitable  for  preclinical  models  of 
breast  cancer  and  clinical  trials 

4.a.  Y1Q1-Y1Q2.  Low  field  MRI:  pre-clinical  development  will  be  done  based  on  our  existing  hardware  that 
was  developed  during  our  preliminary  work.  The  work  will  focus  on  developing  and  construction  of  multiple  RF 
coils  intended  for  optimal  direct  and  indirect  detection  of  hyperpolarized  contrast  agents  in  vivo.  Additional  time 
will  be  spent  on  implementation  of  ultra-fast  MRI  sequences. 

4.b.  Y1Q3-Y1Q4.  Low  field  MRI:  pre-clinical  optimization  will  be  conducted  to  optimize  developed  hardware 
and  MRI  sequences  for  specific  metabolic  contrast  agents  and  their  metabolites  in  vivo. 
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Appendix  2:  Performance  Evaluation  for  QPC  Laser  ES-6507-Z002,  SN  1209201 

Performance  Evaluation  for  QPC  Laser  ES-6507-Z002,  SN  1209201 
(supplement  to  QPC  Product  Test  Report) 


Laura  Walkup,  Brogan  Gust,  Boyd  Goodson 

Department  of  Chemistry  and  Biochemistry ,  Southern  Illinois  University ,  Carbondale,  IL  62901 


Internal  Document  For:  Hyperpolarized  Xenon  Technology  Consortium  (HXTC)  & 

QPC  Lasers 


9  November  2012 


Version  1.0 

(Note:  This  document  will  be  amended  once  the  twin  TEC-cooled  side-arms  and  retro-reflection  become  operational,  and  once 

1”  optics  have  been  tested) 


Reference  Information: 

Nominal  Current 

52.5  A 

Nominal  Max  Current  (Overdrive) 

55  A 

Nominal  Max  Diode  Temperature 

35  °C 

18 


Performance  without  optical  train — 


Current  versus  power  for  various  LDA  temperatures: 


C/D 

cd 

5 


CD 

£ 

O 

CL 


15  20  25  30  35  40  45  50  55 


Current  (amps) 


Water  chiller  temperatures  (°C)  to  achieve  desired  LDA  temp  at  various  currents: 


LDA 

(°C) 

15  amps 

20  amps 

25amps 

30amps 

35amps 

40amps 

45amps 

50  amps 

52. 5  amps 

12 

11.0 

10.2 

9.0 

8.2 

7.8 

6.2 

5.2 

na* 

na* 

18 

17.0 

15.8 

15.0 

14.0 

13.0 

12.0 

11.0 

10.0 

9.3 

20 

18.5 

17.8 

16.8 

16.0 

15.0 

14.0 

13.0 

12.0 

11.5 

22 

20.5 

19.6 

18.6 

18.0 

17.0 

16.0 

15.0 

14.0 

13.5 

25 

23.0 

22.5 

21.5 

20.5 

19.5 

18.5 

17.5 

16.5 

17.1 

28 

26.0 

25.0 

24.0 

23.5 

22.2 

21.2 

20.4 

19.6 

18.9 

30 

28.0 

27.0 

26.0 

25.0 

24.0 

23.3 

22.3 

21.3 

20.8 

*-  lowest  setting  on  water  chiller  is  5°C 


Output  center  wavelength,  as  a  function  of  LDA  temp  and  current: 


LDA 

(°C) 

15  amps 

20  amps 

25amps 

30amps 

35amps 

40  amps 

45amps 

50  amps 

52. 5  amps 

12 

793.18 

793.27 

793.41 

793.54 

793.67 

793.81 

793.94 

na 

na 

18 

793.54 

793.63 

793.74 

793.87 

794.01 

794.14 

794.29 

794.45 

794.53 

20 

793.63 

793.74 

793.87 

793.98 

794.14 

794.27 

794.43 

794.58 

794.61 

22 

793.76 

793.83 

793.98 

793.12 

794.25 

794.39 

794.56 

794.70 

794.74 

25 

793.94 

794.05 

794.16 

794.29 

794.43 

794.56 

794.69 

794.89 

794.91 

28 

794.07 

794.21 

794.32 

794.47 

794.61 

794.74 

794.91 

795.05 

795.00 

30 

794.21 

794.34 

794.45 

794.56 

794.72 

794.87 

795.05 

795.20 

795.27 

[Nominal  resonance:  794.767  nm] 
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notes 


LDA  12  °C 


current 

center 

wavelength 

power 

(W) 

FWHM 

15 

793.18 

10.4 

20 

793.27 

34.8 

25 

793.41 

61.3 

30 

793.54 

88.6 

35 

793.67 

115 

0.27 

40 

793.81 

142 

0.31 

45 

793.94 

169 

0.35 

power  (Watts) 

■  LDA12 

/ 

/' 

( 

)  10  20  30  40  50 

current  (amps) 

LDA  18  °C 


wavelength  (nm) 


notes 


current 

center 

wavelength 

power 

(W) 

FWHM 

15 

793.54 

9.70 

20 

793.63 

33.7 

25 

793.74 

60.0 

30 

793.87 

86.6 

35 

794.01 

113 

0.26 

40 

794.14 

139 

0.31 

45 

794.29 

165 

0.32 

50 

794.45 

189 

0.40 

LDA  20  °C 


notes 


current 

center 

wavelength 

power 

(W) 

FWHM 

15 

793.63 

9.19 

20 

793.74 

33.2 

25 

793.87 

59.3 

30 

793.98 

85.4 

35 

794.14 

111 

0.22 

40 

794.27 

137 

0.23 

45 

794.43 

162 

0.31 

50 

794.58 

187 

0.36 

51 

794.61 

192 

0.37 

LDA  20  °C 


-  15amps 

-  20amps 
25amps 
30amps 
35amps 
40amps 
45amps 

-  50amps 

-  51  amps 
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LDA  22  °C 


notes 

current 

center 

wavelength 

power 

(W) 

FWHM 

15 

793.76 

8.7 

20 

793.85 

32.1 

25 

793.98 

58.1 

30 

794.12 

83.9 

35 

794.25 

110 

0.26 

40 

794.39 

135 

0.27 

45 

794.56 

160 

0.31 

50 

794.70 

185 

0.31 

51 

794.74 

190 

0.37 

52 

794.76 

195 

0.4 

Power  (Watts 

1  — |§i  LDA22 

/ 

■ 

■ 

■ 

/ 

0  10  20  30  40  50 

Current  (amps) 

LDA  25  °C 


notes 


current 

center 

wavelength 

power 

(W) 

FWHM 

15 

793.94 

7.4 

20 

794.05 

30.7 

25 

794.16 

56.3 

30 

794.29 

81.9 

35 

794.43 

107 

0.23 

40 

794.56 

132 

0.23 

45 

794.69 

157 

0.31 

47 

794.78 

167 

0.26 

50 

794.89 

181 

0.36 

— LDA25 

200 


1 

5  ioo- 


0-1 - , - r- 

0  10 


Current  (amps) 
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LDA  28  °C 


notes 


current 

center 

wavelength 

power 

(W) 

FWHM 

15 

794.07 

6.13 

20 

794.21 

28.9 

25 

794.32 

54.1 

30 

794.47 

79.4 

35 

794.61 

104 

0.22 

40 

794.74 

129 

0.26 

45 

794.91 

153 

0.31 

50 

795.05 

175 

0.31 

•S  100- 


/ 


Current  (amps) 


LDA  30  °C 


notes 


current 

center 

wavelength 

power 

(W) 

FWHM 

15 

794.21 

4.71 

20 

794.34 

27.4 

25 

794.45 

51.9 

30 

794.56 

77.2 

35 

794.72 

102 

0.23 

40 

794.87 

126 

0.29 

45 

795.05 

150 

0.29 

50 

795.20 

173 

0.29 

Current  (amps) 


52.5  amps  (nominal  current) 


notes 


LDA 

center 

wavelength 

power 

13.7* 

794.26 

202 

16 

794.40 

200 

18 

794.53 

198 

20 

794.65 

196 

22 

794.76 

193 

24 

794.91 

190 

26 

795.00 

187 

28 

795.16 

184 

22 


wavelength  (nm) 


30 


795.27 


180 


^maximum  chiller  capacity  (LDA  temperature  is  13.7 
°C  at  52.5  amps  when  chiller  is  set  to  5  °C) 
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“Resonance” 


LDA 

current 

power 

center  X 

FWHM 

21.9 

52.5 

194 

794.74 

0.35 

22.8 

51 

188 

794.74 

0.33 

23.1 

50 

182 

794.74 

0.31 

25.7 

45 

157 

794.74 

0.29 

28.4 

40 

132 

794.78 

0.28 

30.9 

35 

99 

794.78 

0.26 

31.7 

33 

91 

794.78 

0.24 

32.2 

32 

84 

794.74 

0.23 

32.8 

31 

78 

794.74 

0.25 

33.1 

30 

68 

794.78 

0.26 

34.1 

28 

62 

794.74 

0.20 
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To  mimic  QPC  specs  but  without  optical  train  mounted  (QPC  specs  include  optical  train) 
Water  chiller  set  at  13°C 

Power  measured  directly  out  of  laser  fibre  (no  optical  train) _ 


LDA  temp  (°C) 

Current  (amps) 

Voltage  (volts) 

Power  (Watts) 

14.1 

15 

13.34 

10.2 

15.0 

20 

13.51 

34.4 

15.9 

25 

13.66 

60.8 

17.0 

30 

13.80 

87.2 

18.1 

35 

13.93 

113 

19.1 

40 

14.06 

133 

20.3 

45 

14.17 

162 

21.2 

50 

14.29 

186 

To  reproduce  QPC  specs  with  2”  optical  train  mounted,  water  chiller  set  to  13°C 
Power  measured  at  “front”  position  of  optical  cell  (-9.85”  from  end  of  optical  train) 


LDA 

temp 

(°C) 

current 

(amps) 

voltage 

(volts) 

QPC  spec 
power 

(W) 

front 

power 

(W) 

dump 

power 

(W) 

power  sum 
(Watts) 

split 

front/dump 

14.1 

15 

13.34 

9.5 

9.48 

0.20 

9.68 

98/2 

15.0 

20 

13.50 

32.0 

32.0 

0.83 

32.83 

97/3 

16.0 

25 

13.65 

57.0 

56.2 

1.58 

57.78 

97/3 

17.1 

30 

13.79 

80.8 

81 

2.34 

83.34 

97/3 

18.1 

35 

13.93 

106.2 

105 

3.14 

108.14 

97/3 

19.1 

40 

14.05 

129.2 

128 

4.06 

132.06 

97/3 

20.1 

45 

14.17 

151.2 

150 

5.09 

155.09 

97/3 

21.2 

50 

14.29 

172.0 

171 

6.00 

177 

97/3 

21.7 

52.5 

14.34 

180 

181 

6.42 

187.42 

97/3 

Power  with  and  without  optical  train  mounted  (to  approximate  power  loss  to  optics) 


W/O  OT 

W/OT 

current 

(amps) 

Power 

(Watts) 

front 

power 

(W) 

dump 

power 

(W) 

power 

sum 

(Watts) 

Difference 

%-age 

power 

lost 

15 

10.2 

9.48 

0.20 

9.68 

0.52 

5.1 

20 

34.4 

32.0 

0.83 

32.83 

1.57 

4.6 

25 

60.8 

56.2 

1.58 

57.78 

3.02 

5.0 

30 

87.2 

81 

2.34 

83.34 

3.86 

4.4 

35 

113 

105 

3.14 

108.14 

4.86 

4.3 

40 

133 

128 

4.06 

132.06 

0.94 

0.7 

45 

162 

150 

5.09 

155.09 

6.91 

4.3 

50 

186 

171 

6.00 

177 

9 

4.8 

52.5 

194 

181 

6.42 

187.42 

6.58 

3.4 
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Total  effective  efficiency 


W/O  OT 

W/OT 

current 

(amps) 

Power 

(Watts) 

front 

power 

(W) 

dump 

power 

(W) 

power 

sum 

(Watts) 

Total  effective  efficiency 
(%-age) 

15 

10.2 

9.48 

0.20 

9.68 

92.9 

20 

34.4 

32.0 

0.83 

32.83 

93.0 

25 

60.8 

56.2 

1.58 

57.78 

92.4 

30 

87.2 

81 

2.34 

83.34 

92.9 

35 

113 

105 

3.14 

108.14 

92.9 

40 

133 

128 

4.06 

132.06 

96.2 

45 

162 

150 

5.09 

155.09 

92.6 

50 

186 

171 

6.00 

177 

91.9 

52.5 

194 

181 

6.42 

187.42 

93.3 
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Beam  diameter  through  optical  train: 


■ - r  p_50 

- - 9.55 - - 


From  XeUS  schematics, 
distances  from  end  of  optical 
train  to  the  OP  cell: 


http  ://www.youtube.com/watch?v=gIE-SPl  Rams 


to  cell  front:  9.85” 
to  cell  middle:  14.68” 
to  cell  back:  19.5” 


Optical  train  mounted  to  laser 
module  in  the  horizontal 
orientation. 

Beam  dumping  in  this 
orientation: 

Retroreflection  | 

Dump  beam  j 

At  52.5  amps  (highest 
operating  current),  dump 
beam  measured  as  6.42  W. 

One  included  heat  sink  was 
mounted  for  retroreflection;  a 
second  heat  sink  has  been 
ordered. 


Beam  collimation  test 
through  2’  ’  optical  train  on 
YouTube 


LDA  16°C,  56 W  forward, 
-1.5W  dump  beam 
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Card  marked  with  2”  diameter  circle: 
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Conclusions 


•  Laser  exhibits  excellent  lineshape  lock  quality  over  nearly  the  entire  range  of  operating  conditions. 

•  Beam  quality  and  collimation  through  optical  train  are  excellent  throughout  relevant  operating  condition  range. 

•  All  required  specs  (in  terms  of  FWHM,  centroid  wavelength  tuning  range,  output  power,  operating  temperature  range,  and 
efficiency  of  optics  train)  are  met. 

•  Heating  of  optical  train  is  not  significant—  even  during  periods  of  prolonged  operation. 

•  Not  yet  tested:  CP  fidelity  of  output;  operation  under  retroreflection;  operation  of  heatsinks  for  cooling  under  retroreflection; 
1”  optics;  long-term  stability  of  output  characteristics. 
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Appendix  3:  Abstracts  Presented 

Abstract  of  Oral  Presentation  #2,  reference  #  (5) 


PHIP  polarizing  and  imaging  strategies  for  13C-lactate  and  lsN-choline 

Eduard  Y.  Chekmenev1,2,3,4,5,  Kevin  W.  Waddell1,2,  Roman  V.  Shchepin1,2,  Aaron  Coffey1,3 

institute  of  Imaging  Science,  and  Departments  of  2Radiology  and  Radiological  Sciences,  3Biomedical  Engineering,  4Biochemistry,  5Vanderbilt-Ingram  Cancer  Center, 

Vanderbilt  University,  Nashville,  TN  37232  United  States 
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Q~  catalyst,  60  °C  ^ 
1 0  atm  para-H2 


203F 


1-13C-phosphoenolpyruvate 

(PEP) 

u  Hyperpolarized 


DL-1  -13C-phospholactate 
(PLAC) 


Introduction:  One  of  the  main  bottlenecks  for  preclinical  and  clinical  application  of  hyperpolarized  MR  is  the  development  and  validation  of  relevant  contrast  agents 
that  can  probe  biochemical  pathways  in  vivo.  Dynamic  nuclear  polarization  (DNP)1  has  been  most  widely  used  to  date  with  the  main  drawback  of  long  (~  1  h) 
polarization  cycles.  Parahydrogen  induced  polarization  (PHIP)  offers  significantly  faster 
preparations  with  hyperpolarization  cycles  as  short  as  1  minute.  However,  it  is  limited  by  the 
availability  of  the  required  unsaturated  molecular  precursors  that  are  necessary  for  molecular 
addition  of  parahydrogen2,3,  which  acts  as  a  source  of  spin  order.  An  additional  requirement  for 
increased  relaxation  time  of  the  hyperpolarized  13C  or  15N  site  is  the  absence  of  directly  attached 
protons.  As  a  result,  the  smallest  PHIP  moiety  consists  of  an  unsaturated  C=C  or  C=C  bond 
adjacent  to  a  labeled  13C  or  15N  site.  For  13C  hyperpolarized  compounds,  this  represents  a  three 
carbon  limitation  successfully  exemplified  by  acrylate  moiety  by  a  number  of  groups.4,5  In 
contrast,  a  leading  DNP  hyperpolarized  metabolic  agent  l-13C-pyruvate  is  also  a  three-carbon 
molecule  carrying  an  extra  oxygen  atom  in  addition  to  the  three-carbon  skeleton  of  propionate. 

Hyperpolarized  l-13C-pyruvate  is  not  amenable  by  PHIP.  Here  we  present  a  strategy  of  using  - 
OH  phosphate  protected  precursor  for  preparation  of  PHIP  hyperpolarized  contrast  agents. 

Nuclear  spin  polarization  of  hyperpolarized  contrast  agents  is  independent  of  the  applied 
magnetic  field  when  detected  in  vivo.  As  a  result,  low  field  MR  detection  of  hyperpolarized 
contrast  agents  in  vivo  can  potentially  be  as  sensitive  as  the  one  at  high  field6  as  long  as  the 
sensitivity  of  detection  RF  coil  is  compensated  with  increased  number  of  turns.  We  recently 
described  0.0475  T  dual  channel  RF  coil  capable  of  sensitive  in  situ  detection  13C  detection7  of 
PHIP  hyperpolarized  contrast  agents.8  First  imaging  results  will  be  presented  using  this  method 
adapted  for  preclinical  molecular  imaging  of  cancer  in  mice. 


7  pmoles 
at  0.0475 

Tesla 

1  scan 

-TTJTTTTJTTT1 

-2  -1 

l"rrl"lll" 

0  1  2 

13C  Frequency  (kHz) 

0  10  20  30  40 

Time  (seconds) 


Methods:  We  used  PHIP  polarizer  with  in  situ  detection  capability  previously  described8  with  the 
exception  that  we  utilized  a  new  dual  channel  RF  coil  optimized  for  13C  PHIP 
detection.  97%  parahydrogen  was  produced  using  in-house  built  parahydrogen 
generator.9  Catalyst  and  PHIP  precursor  preparations  were  done  according  to 
previously  described  protocols.  l-13C-phosphoenolpyruvate  (Isotec)  was  used  as 
PHIP  precursor  for  hyperpolarization  of  l-13C-phospholactate,  Fig.  1. 


Results:  PHIP  polarization  P  ~  1%  was  achieved  with  l-C-phospholactate,  Fig. 
1.  13C  77  of  36  s  was  measured  in  situ  at  0.0475  T,10  which  is  suitable  for  in  vivo 
application,  Fig.  1. 

The  NMR  probe  head  for  XH  and  X  nucleus  (X  =  13C  or  15N)  was  constructed  for 
operation  at  a  static  field  of  B0  =  0.0475  T.  The  RF  coil,  Fig.  2a,  consists  of  two 
single  channel  circuits  sharing  a  ground  with  capacitors  for  achieving  suitable 
tuning  range  and  impedance  matching  for  each  channel.  The  NMR  probe 
frequencies  of  interest  are  2.02  MHz  (1H),  0.508  MHz  (13C),  and  0.205  MHz  (15N) 
at  the  Halbach  array's  0.0475  T  B0  field  strength.  The  1 H  circuit  consists  of  an 
outer  16  turn  (8  turns  each  side)  70  mm  x  130  mm  (I.D.  x  length)  Helmholtz 
saddle  coil.  The  second  X  channel  coil  is  a  206  turn  50  mm  x  170  mm  (I.D.  x 
length)  single  layer  solenoid  closely  fitted  to  the  high-pressure  reactor  used  for 
molecular  addition  of  parahydrogen  to  unsaturated  molecular  precursors.  HEP  (2- 
hydroxyethyl  propionate- 1-13C, 2, 3, 3-d3)  was  polarized  with  PASADENA3  by 
catalytic  hydrogenation  from  HEA  (2-hydroxyethyl  acrylate-l-13C,2,3,3-d3)  using 
up  to  97  %  enriched  parahydrogen.  Fig.  2b  shows  in  situ  13C  detection  of  15 
micromoles  (<  2  mg)  of  HEP  with  20%  polarization,  or  signal  and  polarization 
enhancement  s  ~  5,000,000,  with  a  13C  SNR  ~  2,  300  and  FWHM  =  26  Hz. 
Furthermore,  the  13C  MR  sensitivity  of  the  presented  0.0475  T  RF  probe  was 


Figure  1.  a)  PHIP  polarization  of  l-C-PEP  resulting  in  phosphate 
protected  l-13C-lactate,  l-13C-phospholactate,  b)  single  scan  13C  spectrum 
of  7  micromoles  of  1%  hyperpolarized  l-13C-phospholactate  recorded  in 
situ  of  0.0475  T  PHIP  polarizer,  c)  13C  T,  decay  of  hyperpolarized  13C- 
phospholactate  at  0.0475  T. _ 


°  ^  Boz  ~ 

'sj  Wm*) 

^^I^^MX-COil)  X  < 

Magnet  Bore  * 


b 

^.150' 

\  100' 
J 

£=  50  • 
0  ■ 


*200^ 

15  micromoles  13C 
1  scan 

4.IO-6 
=  4.7  T 


504  506  508  510  512 

13C  Larmor  Frequency  (kHz) 


4.IO-6 

B0  =  47.5  mT 


-2  0  2  4 

13C  Chemical  Shift  (kHz) 


-4-2  0  2 

13C  Chemical  Shift  (kHz) 


Figure  1.  a)  Alignment  of  *H-X  RF  coils,  magnetic  field  and  magnet  bore,  b)  Single 
acquisition  spectrum  of  15  micromoles  (<  2  mg)  of  hyperpolarized  HEP  contrast  agent 
with  20%  polarization  or  enhancement  s  ~  5,000,000  at  0.0475  T,  c)  &  d)  sensitivity 
comparison  at  high  and  low  field.  13C  spectroscopy  of  1.0  g  of  sodium  l-13C-acetate 
solution  in  2.8  mL  99.8%  D20. 


compared  to  that  of  a  commercial  4.7  T  small-animal  RF  coil  with  smaller  volume.  We  used  a  sample  of  sodium  1-  C-acetate  (1  g  in  2.8  mL  D20  solution)  polarized  to 
P  ~  4*10'6  for  single  scan  13C  acquisition,  Fig.  2  using  spectra  width  of  10  kHz.  The  high  field  (4.7  T)  volume  coil  yielded  an  SNR  of  120  with  a  FWHM  of  6  Hz,  while 
the  13C  solenoid  coil  of  the  0.0475  T  H-X  probe  achieved  an  SNR  of  28  with  a  FWHM  of  25  Hz. 


Discussion:  The  method  of  -OH  phosphate  protected  and  stabilized  C=C  bond  offers  an  opportunity  to  polarize  metabolic  contrast  agents  not  amenable  in  the  past. 
Presented  here  is  the  preliminary  result  of  PHIP  hyperpolarization  of  l-13C-phospholactate.  PHIP  hydrogenation  is  efficient  and  blood  phosphatases  cleave  phosphate 
moiety  after  tracer  injection  in  living  organism  (data  not  shown)  producing  metabolically  relevant  (l-13C-lactate)  contrast  agent.  The  method  can  be  potentially  extended 
to  protected  15N-cholines.u  13C  hyperpolarized  lactate  is  a  potent  metabolic  imaging  contrast  agent.12 

Conclusion:  MR  sensitivity  compensation  with  multi-turn  RF  coils6,8  offers  an  opportunity  for  sensitive  MRI  of  hyperpolarized  contrast  agents  in  low  magnetic  fields, 
where  magnetic  susceptibility  and  subject  induced  Bj  inhomogeneity  can  be  significantly  minimized  to  the  extent  that  no  scan  preparation  (e.g.  shimming  and  B, 
mapping)  is  required  potentially  allowing  to  true  sub-second  scan  speed  of  hyperpolarized  MR  with  minimized  cost  of  MRI  scanner. 
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Owing  to  the  detection  sensitivity  provided  by  their  high,  non-equilibrium 
magnetization,  hyperpolarized  (HP)  gases  (e.g.  129Xe  and  3He)  are  finding  use  in  a 
growing  number  of  MRS/MRI  applications — ranging  from  biomedical  imaging  and 
spectroscopy,  to  probing  molecular  and  materials  surfaces  [1-4].  Although  3He  has 
a  nearly  three-fold  greater  gyromagnetic  ratio  (and  until  recently  [5-8],  was 
generally  considered  easier  to  hyperpolarize),  129Xe  is  more  naturally  abundant  (by 
several  orders  of  magnitude),  and  xenon  possesses  higher  solubility  in  blood  and 
other  tissues.  Moreover,  129Xe’s  proclivity  for  interacting  with  other  substances  and 
its  much  wider  chemical  shift  range  make  it  a  more  sensitive  MR  probe  of  chemical 
and  biological  environments  [1]. 

HP129Xe  is  usually  created  by  spin-exchange  optical  pumping  (SEOP)  [9].  It 
is  traditionally  expected  that  high  129Xe  polarizations  (PXe)  can  only  be  obtained  with 
low  in-cell  Xe  densities,  because:  (1)  higher  Xe  densities  increase  the  destruction  of 
the  alkali  metal  polarization  from  non-spin-conserving  collisions;  and  (2)  higher  total 
pressures  tend  to  quench  the  (more  efficient)  van  der  Waals  contribution  to  Rb-Xe 
spin  exchange  [9].  Indeed,  many  polarizer  designs  tend  to  go  to  great  lengths  to  produce  large  amounts  of  HP129Xe  while 
still  satisfying  this  condition  within  the  cell.  In  past  work  from  the  SlUC/Nottingham  collaboration  [10,11]  exploring  batch¬ 
mode  and  “stopped-flow”  [12]  Rb/Xe  SEOP  under  conditions  of  high  resonant  laser  flux,  an  unexpected,  inverse 
relationship  was  found  between  the  optimal  temperature  for  performing  Rb/Xe  SEOP  and  the  in-cell  Xe  density  that  could 
be  exploited  to  provide  surprisingly  high  PXe  values  at  high  in-cell  Xe  densities  [10,11]. 

Scaling  up  from  these  results,  our  consortium  recently  constructed  an 
automated  large-scale  (>1  L/hr)  first-generation  ‘open-source’  129Xe  polariser 
(dubbed  “XeNA”).  Comprised  mostly  of  off-the-shelf  components  (including  a 
200  W  VHG-narrowed  LDA  laser),  XeNA  can  be  used  for  clinical,  pre-clinical, 
and  materials  NMR/MRI  applications — and  is  now  cleared  by  FDA  for 
operation  at  Brigham  &  Women’s  Hospital.  Unlike  most  clinical-scale  Xe 
polarizers,  this  first-generation  device  runs  with  Xe-rich  gas  mixtures  (up  to 
2000  Torr  in  a  500  cc  cell)  in  either  single-batch  or  stopped-flow  mode,  in  part 
negating  the  usual  requirement  to  cryo-collect  the  HP129Xe.  In-cell  (during 
SEOP)  PXe  values  of  up  to  -90%,  -57%,  -50%,  &  -28%  have  been  measured 
for  Xe  partial  pressures  of  -300,  -500,  -760,  &  -1570  Torr,  respectively;  PXe 
of  -41%  and  -28%  have  been  measured  after  transfer  of  the  HP129Xe  to 
Tedlar  bags  at  3  T  [13]. 

Taking  what  we  have  learned  from  constructing  our  first-generation  129Xe  polarizer,  we  have  designed  (and  are 
currently  constructing)  our  second-generation  129Xe  polarizer  (dubbed  “XeUS”).  XeUS’s  design  includes  a  variety  of 
improvements  and  new  technologies.  For  example,  the  next-generation  water-cooled  200  W  VHG-narrowed  LDA  laser 
possesses  a  novel  attached  air-cooled  optical  train  assembly  (designed  by  QPC  and  the  consortium).  Implementation  of  a 
large  thermoelectric  cooler  (TEC)  and  quiet  miniature  air  compressor  eliminates  the  need  for  external  gas  or  liquid-N2 
supplies  for  heating  /  cooling  of  the  OP  cell  oven  or  for  pneumatic  valve  operation;  the  polycarbonate  oven  for  the  OP  cell 
is  manufactured  via  3D  printing  and  houses  the  TEC  heating  and  cooling  of  the  cell.  The  GUI-controlled  and  automated 
gas  manifold  and  its  operation  is  simplified  in  part  by  the  use  a  single  tank  of  pre-mixed  129Xe/N2  gas  and  maintaining 
overpressure  throughout  the  space  upstream  from  the  cell  (a  change  that  should  also  prolong  the  operational  lifetime  of 
the  alkali  metal  in  the  OP  cells  by  reducing  oxidation,  while  also  reducing  the  delay  time  between  the  preparation  of 
different  HPXe  batches).  The  four  magnetic-coil  configuration  follows  a  “Barker”  design  [14]  providing  a  more 
homogeneous  magnetic  field  compared  to  a  regular  Helmholtz  setup.  The  device  is  expected  to  produce  >1  L/hr  of  129Xe 
with  PXe  >20%.  Importantly,  XeUS  will  be  integrated  with  an  apparatus  that  will  utilize  the  resulting  HPXe  as  a 
hyperpolarization  source  for  the  nuclear  spins  of  other  species  (“Xenon-Induced  Polarization”,  or  XIP). 

Acknowledgments:  Work  at  SIUC  and  Vanderbilt  is  supported  by  a  DoD  CDMRP  Era  of  Hope  Award  W81 XWH-1 2-1 -01 59/BC1 12431.  M.J.B  is  supported  by  the 
School  of  Medical  &  Surgical  Sciences,  U.  of  Nottingham.  References:  [1]  Goodson,  JMR,  155  157  (2002);  [2]  Mugler,  et  al.,  Magn.  Reson.Med.  37  809  (1997),  [3] 
Leawoods,  et  al.,  Concepts  Magn.  Reson.  13  Til  (2001);  [4]  Patz,  et  al.,  EurJ.  Radiol.  64  335  (2007),  [5]  Driehuys,  et  al.  APL  69, 12  (1996);  [6]  Ruth,  et  al.  APB  68, 
93  (1999);  [7]  Zook,  et  al.,  JMR  159,  175  (2002);  [8]  Ruset  et  al.,  PRL  96,  053002  (2006);  [9]  Walker  &  Happer,  Rev.  Mod.  Phys.  69,  629  (1997);  [10]  Whiting, 
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Owing  to  the  detection  sensitivity  provided  by  their  high,  non-equilibrium  nuclear  spin  polarizations,  hyperpolarized  (HP) 
noble  gases  (e.g.  129Xe  and  3He)  are  utilized  in  a  growing  number  of  MRS/MRI  applications — ranging  from  biomedical 
imaging  and  spectroscopy  to  probing  molecular  and  materials  surfaces  [1].  Lung  imaging  with  HP129Xe  is  of  particular 
interest  [2];  although  the  gyromagnetic  ratio  of  129Xe  is  smaller  than  that  of  3He,  xenon  is  far  more  naturally  abundant  (and 
possesses  higher  solubility  in  blood).  However,  a  major  obstacle  towards  implementing  HP129Xe  for  clinical  imaging  has 
been  the  ability  to  cheaply  and  easily  generate  large  quantities  of  HP  gas  with  high  129Xe  polarizations  (PXe).  HP129Xe  is 

usually  created  via  spin-exchange  optical  pumping  (SEOP  [3]).  It  is  generally 
anticipated  that  high  PXe  values  are  achievable  via  SEOP  only  with  low  in-cell 
Xe  densities;  thus,  most  large-scale  xenon  polarizing  devices  (“polarizers”, 
including  those  from  commercial  sources)  go  to  great  lengths  to  satisfy  this 
condition  within  the  cell.  Such  polarizers  also  tend  to  be  complex,  expensive, 
and  strictly  regulated  by  their  controlling  entities — factors  that  hinder  large 
scale  HP^9Xe  generation  for  many  unaffiliated  laboratories. 

Building  upon  previous  results  exploring  129Xe  SEOP  under  conditions  of  high 
resonant  laser  fluxes  and  high  xenon  densities  [4,5],  here  we  present  results 
obtained  with  our  collaboration’s  first-generation  large-scale  (>1  L/hr)  “open- 
source”  xenon  polarizer  for  clinical,  pre-clinical,  and  materials  NMR/MRI 
applications.  Comprised  of  mostly  off-the-shelf  components  (including  a  200 
W  VHG-narrowed  LDA  laser),  our  automated,  modular  polarizer  is  easy  to 
use  and  employs  a  simplified  “open-source”  design,  readily  implementable  in 
other  laboratories.  Unlike  most  clinical-scale  Xe  polarizers,  our  device  runs 
with  xenon-rich  gas  mixtures  (up  to  1800  Torr  in  a  500  cc  cell)  in  either  single¬ 
batch  or  stopped-flow  mode,  in  part  negating  the  usual  requirement  to  cryo- 
collect  the  HP129Xe  (and  potentially  opening  a  door  to  studies  of  quadrupolar 
noble  gas  nuclei  [6]).  Up  to  four  independent  methods  were  utilized  to  probe 
the  spin  polarization  achieved  with  the  polarizer — including  in  situ  field-cycled 
near-IR  spectroscopy  [4]  (to  probe  the  alkali  metal  electron  spin  polarization), 
in  situ  low-field  129Xe  NMR  (calibrated  with  thermal  NMR);  xenon  transfer 
to  47.5  mT  for  129Xe  NMR  (calibrated  with  thermal  13C  NMR),  and  xenon 
transfer  and  129Xe  NMR  at  3  T  (calibrated  against  thermal  129Xe  NMR)— and 


o  400  eoo  1200  i&ao 
Xe  Partial  Pressure  (Torr) 

Fig.  1  {lefty  Pxe  values  measured  at  5.26  mT, 
47.5  mT,  and/or  3  T,  plotted  versus  xenon 
partial  pressure  during  SEOP.  Labels  'b.t.’  and 
'a.t.’  respectively  refer  to  measurement 
obtained  ’’before”  and  ’’after”  Xe  gas  was 
transferred  to  another  container.  Inset:  2D 
FLASH  129Xe  projection  image  of  HPXe  in  a 
Tedlar  bag  following  transport  to  3  T  (FOV:  40 
cm;  slice  thickness:  1.7  cm;  TE/TR:  2.6/5. 1  ms, 
flip  angle:  ~4°;  SNR  -40:1). 


excellent  agreement  was  observed  among  these  different  approaches.  In-cell 
(during  SEOP)  PXe  values  of  up  to  ~90%,  ~57%,  -50%,  and  ~25%  have  been  measured  for  Xe  partial  pressures  of  ~300, 
~500,  -760,  and  -1570  Torr,  respectively  (Fig.  1).  PXe  values  of  -41%  and  -28%  have  been  measured  after  transfer  of  the 
HP129Xe  to  Tedlar  bags  at  3  T  (Fig.  1,  inset).  In-bag  TVs  of  -38  min  and  nearly  6  hr  at  low  field  and  3  T,  respectively, 
have  been  recorded — more  than  long  enough  for  a  variety  of  potential  experimental  modalities.  The  device  has  received 
FDA  and  internal  IRB  approval.  Finally,  we  will  report  on  in  vivo  experiments,  which  are  now  underway. 
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Hyperpolarization  technologies  have  lead  to  sensitivities  greater  than  conventional  MRI  or  MRS  by  several  orders  of 
magnitude  increasingly  improving  ability  to  assay  cellular  metabolism  at  the  molecular  level1,2.  Clinical  implementation  has 
hinged  on  two  critical  developments:  on-site  production  capacity  for  short-lived  contrast  agents  and  capability  for  efficient 
MRI  of  these  agents3.  While  the  increasing  commercial  availability  of  DNP,  SABRE,  PASADENA,  and  PHIP  based 
systems4  have  made  contrast  agents  more  accessible,  the  latter  remains  largely  unaddressed  in  the  low  field  regime 
outside  the  domain  of  hyperpolarized  noble  gas  imaging.  The  success  of  low  field  noble  gas  imaging  systems5,6  has 
demonstrated  the  utility  of  detection  at  low  magnetic  field  strengths.  Hyperpolarized  technologies  permit  preparation  of 
nuclear  spin  states  where  the  detected  nuclear  spin  polarization  is  independent  of  the  applied  magnetic  field  B0  and  signal 
levels  depend  on  the  magnetic  moment  of  the  observed  nucleus  and  Larmor  frequency  co0 ■  This  NMR  detected  signal 
voltage,  known7,8  to  correspond  to  the  electromotive  force  ( Emf)  arising  from  Faraday’s  law  of  induction,  is  linearly 
proportional  to  the  number  of  turns  (A/)  of  RF  coil  and  the  change  of  magnetic  flux  Oe  through  a  single  loop.  The  Emf  may 
also  be  expressed  as  \Emf\  N-  co0  •  P  •/!,  where  Pis  the  nuclear  spin  polarization,  and  g  is  the  magnetic  moment  of 


nuclear  spin9.  Multi-turn  coils  at  low  co0  creating  high  levels  of  induced  Emf  can  provide  a  highly  sensitive  means  for  MR 
detection  of  hyperpolarized  compounds.  The  practicality  of  this  approach  to  achieving  high  sensitivity  is  demonstrated 
here  for  13C  and  ^  imaging  of  sodium  1-13C-acetate  phantoms  with  constant  polarization  simulating  the  condition  of  the 


hyperpolarized  state  at  2.02  MHz  and  0.508  MHz  Larmor  frequencies. 

Two  spherical  phantoms  of  sodium  1-13C-acetate  were  used  for  'H  and 
performed  on  1.0  g  of  sodium  1-13C-acetate  in  2.8  mL  of  99.8%  D20,  with 
acetate  in  16.3  mL  of  99.8%  D20  (Fig.  1 ).  1H  (flip  angle  a 
=  18°,  spectral  width  (SW)  of  10  kHz,  6.4  ms  acquisition 
time)  and  13C(a  =  90°,  SW  =  5  kHz,  6.4  ms  acquisition 
time)  images  were  acquired  on  a  4.7  T  Varian  scanner 
using  a  multinuclear  Doty  Scientific  RF  probe  with  a 
FSSFP  sequence.  The  phantoms  were  similarly  imaged 
on  a  0.0475  T  Magritek  spectrometer  with  a  Magritek 
gradient  coil  insert  and  custom  in-house  multinuclear  RF 
probes  of  same  38  mm  diameter  and  very  similar  coil 
volume  using  a  fast  gradient  echo  sequence  with  the 
same  imaging  parameters  except  for  13C  flip  angle  a  = 

18°,  not  90°.  For  13C,  the  phantom  was  imaged  5  s  after 
pre-polarization  at  7  T,  after  which  13C  polarization  was 
116.0  ±  0.5  %  4.7  T  thermal  equilibrium  polarization.  ^ 
imaging  was  performed  5  s  after  pre-polarization  at  9.4  T 
with  65  ±  1  %  4.7  T  polarization. 

To  the  best  of  our  knowledge,  this  is  the  first 
demonstration  that  low  field  imaging  can  approach  the 
sensitivity  of  high  field  one  under  conditions  of  field 
independent  polarization.  Based  on  these  results,  we 
anticipate  that  low  field  in  vivo  hyperpolarized  MRI  would 
outperform  the  sensitivity  of  the  high  field,  because  of  (i) 
much  lower  magnetic  field  susceptibility  resulting  in  much 
better  B0  field  homogeneity  across  the  subject,  (ii) 
potential  use  of  cryo-cooled  RF,  and  (iii)  more 
homogeneous  B1  and  significantly  reduced  specific 
absorption  rate  (SAR). 
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Figure  1.  2D  images  acquired  on  a  4.7  T  Varian  scanner  with  a 
multinuclear  Doty  RF  coil  and  on  a  0.0475  T  Magritek  spectrometer 
equipped  with  a  gradient  insert  and  in-house  multinuclear  RF  coils 
for  identical  samples  with  identical  imaging  parameters  except  13C 
flip  angle.  13C  imaging  of  5.18  g  of  sodium  1-13C-acetate  in  99.8% 
D20  (A)  at  4.7  T  and  (B)  at  0.0475  T  5  s  after  pre-polarization  at  7 
T.  imaging  of  1.0  g  of  sodium  1-13C-acetate  in  99.8%  D20  (C)  at 
4.7  T  and  (D)  at  0.0475  T  5  s  after  pre-polarization  at  9.4  T. 
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Hyperpolarization  of  13C  Tracers  for  Molecular  Imaging  of  Cancer 

Panayiotis  Nikolaou,  Eduard  Chekmenev 


Background 

Our  collaboration  (Rosen  -  Harvard  University,  Goodson  -  Southern  Illinois  University,  Patz  -  Brigham  and 
Women’s,  Barlow  -  University  of  Nottingham,  UK)  has  recently  constructed  a  large-scale  (>1  L/hr)  129Xe 
polarizer  (XENA)  for  clinical,  pre-clinical,  and  materials  NMR/MRI  applications  (Fig.  1).  Comprised  mostly 
of  off-the-shelf  components,  our  automated,  modular  polarizer  is  easy  to  use  and  employs  a  simplified 
“open-source”  design  that  will  be  readily  implementable  in  other  laboratories.  Furthermore,  the  2nd 
generation  polarizer  funded  by  the  Department  of  Defense  will  be  constructed  at  VUIIS  and  will  be 
implemented  in  Xenon  Induced  Polarization  experiments  of  13C  tracers  for  molecular  imaging  of  cancer. 


usually  prepared  by  spin-exchange  optical  pumping 
(SEOP)  [1].  It  is  typically  expected  that  high  129Xe  polarizations  (PXe) 
can  only  be  obtained  with  low  in-cell  Xe  densities,  because:  (1)  higher 
Xe  densities  increase  the  destruction  of  the  alkali  metal  polarization 
from  non-spin-conserving  collisions;  and  (2)  higher  total  pressures  tend 
to  quench  the  (more  efficient)  van  der  Waals  contribution  to  Rb-Xe  spin 
exchange  [1].  Indeed,  many  polarizer  designs  produce  large  amounts 
of  HP129Xe  while  still  satisfying  this  low-density  condition  within  the  cell. 
In  recent  work  from  the  SlUC/Nottingham  collaboration  [2,3]  exploring 
batch-mode  and  “stopped-flow”  [3]  Rb/Xe  SEOP  under  conditions  of 
high  resonant  laser  flux,  an  unexpected,  inverse  relationship  was  found 
between  the  optimal  temperature  for  performing  Rb/Xe  SEOP  and  the 
in-cell  Xe  density  [2,3] — an  effect  that  was  recently  exploited  to 
achieve  surprisingly  high  PXe  values  at  high  in-cell  Xe. 

Results 

129Xe  polarizations  of  -70%,  -44%,  and  -30%  (calibrated  with  built-in 
low-field  NMR)  have  been  measured  for  Xe  partial  pressures  of 
300,  760,  and  1600  Torr,  respectively  using  XENA.  XENA  was 
installed  in  Brigham  and  Women’s  Hospital  for  on-going  clinical  trial. 

Significance 

Unlike  most  other  clinical-scale  Xe  polarizers,  this  device  runs  with 
xenon-rich  gas  mixtures  (up  to  2000  Torr  in  a  500  cc  cell)  in  either 
single-batch  or  stopped-flow  mode,  in  part  negating  the  usual 
requirement  to  cryo-collect  the  hyperpolarized  129Xe  (simplifying  gas 
handling  while  eliminating  polarization  losses  associated  with  freeze- 
thaw  cycles).  One  goal  with  the  open-source  design  is  to  facilitate 
greater  implementation  of  HPXe  MR  technology  in  other  labs  and 
feedback  for  future  improvements;  moreover,  shared  designs  in 
multiple  locations  should  also  facilitate  comparisons  of  results,  data 
sharing,  accelerated  development  and  clinical  use.  Our  design  will 
also  allow  exploration  of  alternate  SEOP  approaches  including  optical 
excitation  geometries,  HPXe  storage  and  delivery  methodologies,  and 
the  use  of  cesium. 
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Figure  1:  (top)  Photograph  of  1st 
generation  automated  Xe 
polarizer;  (bottom)  low  field  MR 
spectroscopy  of  produced 
hyperpolarized  129Xe  at  natural 
abundance  of  26.4%  at  300  torr 
Xe  in  external  portable  phantom 
measured  at  47.5  mT  and 
(inset)  in  situ  quality  assurance 
measurement  in  optical 
polarizing  cell  at  5.3  mT  before 
transfer  to  47.5  mT  system 
(almost  no  %Pxe  loss). 
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Appendix  4:  Manuscripts  Published  and  Accepted 

One  peer-reviewed  manuscript  supported  by  this  funding  was  accepted  (7)  and  one  manuscript  was 
conditionally  accepted  by  PNAS(1).  The  copies  of  published  and  accepted  manuscript  are  attached  in  the  end 
of  this  report. 
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Abstract:  Reliable  early  assessment  of  breast  cancer  response  to  neoadjuvant  therapy  (NAT) 
would  provide  considerable  benefit  to  patient  care  and  ongoing  research  efforts,  and  demand  for 
accurate  and  noninvasive  early-response  biomarkers  is  likely  to  increase.  Response  assessment 
techniques  derived  from  quantitative  magnetic  resonance  imaging  (M  Rl)  hold  great  potential 
for  integration  into  treatment  algorithms  and  clinical  trials.  Quantitative  MRI  techniques 
already  available  for  assessing  breast  cancer  response  to  neoadjuvant  therapy  include  lesion 
size  measurement,  dynamic  contrast-enhanced  MRI,  diffusion-weighted  MRI,  and  proton 
magnetic  resonance  spectroscopy.  Emerging  yet  promising  techniques  include  magnetization 
transfer  MRI,  chemical  exchange  saturation  transfer  MRI,  magnetic  resonance  elastography, 
and  hyperpolarized  M  R.  Translating  and  incorporating  these  techniques  into  the  clinical  setting 
will  require  close  attention  to  statistical  validation  methods,  standardization  and  reproducibility 
of  technique,  and  scanning  protocol  design. 

Keywords:  treatment  response,  presurgical  treatment,  neoadjuvant  chemotherapy 


Introduction 

C  hemotherapy  and  hormonal  therapy  for  early  stage  breast  cancer  can  be  admi  nistered 
in  either  the  adjuvant  (after  surgery)  or  neoadjuvant  (before  surgery)  settings.  Potential 
advantages  to  neoadjuvant  therapy  (NAT)  include  presurgical  reduction  of  tumor  bur¬ 
den,  which  may  allow  certain  patients  to  undergo  breast  conservation  therapy  rather 
than  mastectomy,  and  earlier  treatment  of  possibleoccultmicrometastatic  diseasewith 
the  primary  breast  mass  acting  as  a  "marker"  for  treatment  effectiveness  (as  opposed 
to  adjuvant  chemotherapy,  where  no  such  marker  exists  for  response  of  systemic 
micrometastases).1  NAT  also  allows  patients  time  to  undergo  genetic  testing  if  there  is 
a  suspicion  of  an  underlying  BRCA  1/2  mutation,  which  if  found  may  prompt  patients 
to  consider  mastectomy  rather  than  lumpectomy.  At  present,  NAT  is  offered  primar¬ 
ily  to  patients  with  larger  tumors,  tumors  fixed  to  the  chest  wall,  or  clinically  matted 
lymph  nodes  or  skin  involvement.  However,  current  clinical  trials  are  evaluating  the 
use  of  different  NAT  regimens  in  multi  pie  patient  groups,  including  those  with  smaller 
tumors,2  and  it  is  anticipated  that  these  studies  will  lead  to  increased  use  of  NAT  in 
early  stage  breast  cancer  across  a  broad  spectrum  of  patients. 

With  more  breast  cancer  patient  receiving  NAT,  more  available  NAT  regimens, 
and  multiple  emerging  drug  agents  and  combinations  requiring  evaluation  and  com¬ 
parison  with  existing  options,  reliable  assessment  of  treatment  response  has  emerged 
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as  an  important  challenge  in  both  the  clinical  and  research 
environments.  Pathologic  response  has  been  established  as 
an  independent  prognostic  marker  for  overall  survival  in 
breast  cancer  and  i  s  currently  the  gold  standard  for  assessi  ng 
response  to  NAT,3  but  techniques  that  can  provide  reliable 
response  assessment  earlier  in  the  course  of  therapy  are  in 
high  demand,  driven  by  several  considerations.  First,  from 
a  clinical  perspective,  accurate  early  response  assessment 
would  providetheopportunity  to  replacean  i  neffecti  ve  treat¬ 
ment  with  an  alternative  regimen,  and  in  so  doing  potentially 
avoid  or  curtail  debilitating  side  effects  or  toxicities,  such 
as  cardiotoxicity  from  anthracycline  agents  or  neuropathy 
from  taxanes.  Second,  patients  with  disease  that  is  refrac¬ 
tory  to  multiple  NAT  regimens  could  be  referred  directly  to 
surgery,  assuming  this  determination  is  made  early  enough 
that  their  disease  is  still  surgically  resectable.  Third,  from 
a  research  perspective,  accurate  early  response  assessment 
would  allow  for  determination  of  treatment  efficacy  on  a 
much  shortened  timescale,  with  important  ramifications  for 
clinical  trial  design. 

At  present,  there  is  no  uniform  approach  to  the  early 
assessment  of  breast  cancer  response  to  NAT.  Palpation, 
probably  the  most  widely  used  technique  in  the  clinical  set¬ 
ting,  is  inaccurate  for  predicting  pathologic  response45  and  is 
poorly  suited  for  assessment  of  small  tumors.  Serial  biopsy  is 
invasive  and  is  associated  with  sampling  problems  in  hetero¬ 
geneous  tumors.  N  oni  nvasi  ve  i  magi  ng  techni  ques,  including 
mammography,  ultrasound,  and  conventional  M  R I ,  areoften 
deployed  in  the  clinical  setting  as  ad  hoc  problem-solving 
tools  and  in  the  research  setting  as  secondary  endpoints, 
but  no  imaging-based,  early-response  biomarker  has  been 
suitably  validated  as  sufficiently  predictive  of  long-term 
outcomes  to  become  incorporated  either  as  standard  of  care 
in  the  clinical  setting  or  as  a  routine  component  of  all  clini¬ 
cal  trials.  The  stage  is  now  set  for  development  of  accurate, 
noninvasive,  early-response  biomarkers  for  integration  into 
both  the  clinical  and  research  environments. 

This  contribution  discusses  quantitative  magnetic  reso¬ 
nance  imaging  (M  R I )  as  a  promising  platform  from  which 
to  develop  and  deploy  these  biomarkers.  We  begin  by 
discussing  the  motivations  for  using  quantitative  M  Rl  for 
assessing  breast  cancer  response  to  NAT.  We  then  describe, 
with  illustrative  examples,  several  currently  available  quan¬ 
titative  M  Rl  methods  including  lesion  size  measurement, 
dynamic  contrast-enhanced  MRI  (DC E-M  Rl ),  diffusion- 
weighted  M  Rl  (DW-M  Rl),  and  proton  magnetic  resonance 
spectroscopy  (M  RS).  Looking  a  bit  further  on  the  horizon, 
we  then  discuss  magnetization  transfer  (M  T)  M  R  I,  chemical 


exchange  saturation  transfer  (CE  ST)  M  Rl,  M  R  elastography, 
and  hyperpolarized  M  R,  methods  that  have  been  deployed  in 
other  disease  sites  and  that  may  be  particularly  well  suited  to 
breast  cancer  response  assessment.  We  cond  ude  by  address- 
ing  certain  practical  challenges  in  the  clinical  translation  of 
quantitative  M  Rl  methods. 

Rationale  behind  quantitative  MRI 
methods  for  response  assessment 

Cancer  imaging  is  undergoing  a  paradigm  shift  in  which 
quantitative  answers  are  increasingly  being  sought  for  ques¬ 
tions  that  have  historically  motivated  a  qualitative  response. 
The  query  "is  this  cancer  responding  to  therapy"  has  tradi¬ 
tionally  been  answered  by  clinical  imaging  based  on  the 
subjective  impressions  of  the  observer.  However,  modern 
cancer  research  and  treatment  now  requires  objective  and 
reproducible  response  assessment  variables  with  which  to 
evaluate  and  compare  different  treatment  strategies.  For 
evaluating  treatment  response,  qualitative  imaging  interpre¬ 
tation  is  yielding  ground  to  quantitative  imaging  response 
parameters  that  can  be  integrated  with  other  quantitative 
clinical  datasets  for  rigorous  statistical  evaluations.67 

The  current  mainstay  of  objective  imaging-based 
response  assessment  for  solid  malignancies  is  the  Response 
Evaluation  Criteria  in  Solid  Tumors  (RECIST),  a  set  of 
published  guidelines  for  image  acquisition,  lesion  evaluation 
and  measurement,  and  response  categorization.8  R  E  C I  ST  has 
been  successful  in  providing  a  standardized  approach  for 
imaging-based  response  assessment  and  facilitating  "apples 
to  apples"  comparisons  of  different  cancer  treatments,  but 
its  emphasis  on  changes  in  lesion  size  has  been  criticized  as 
failing  to  capture  meaningful  changes  in  tumor  biology.910  In 
particular,  RECIST  may  underestimate  the  antitumor  efficacy 
of  newer  drug  agents  with  cytostatic  rather  than  cytotoxic 
effect,  where  tumor  shrinkage  may  be  minimal  or  quite 
delayed.  The  future  of  imaging-based  quantitative  response 
assessment  will  likely  involve  supplementing  traditional 
size-based  techniques  with  more  advanced  morphological, 
functional,  physiological,  cellular,  and  molecular  imaging 
methods. 

For  a  variety  of  reasons,  MRI  is  a  very  promising  plat¬ 
form  from  which  to  develop  advanced  quantitative  imaging 
biomarkers  for  assessing  breast  cancer  response  to  NAT. 
First,  M  Rl  is  already  well  established  as  a  clinically  useful 
technique  in  breast  cancer  detection  and  characterization,11 
and  its  use  for  response  assessment  would  be  a  natural 
extension  of  its  current  role  in  clinical  care.  Second,  as  an 
intrinsically  digital  technique,  M  Rl  is  capable  of  generating 
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quantitative  datasets  for  direct  entry  into  statistical  analyses, 
thus  bypassing  the  need  for  analog-to-digital  conversion  or 
subjective  user  interpretation  that  can  add  to  measurement 
variability.  Third,  MRI  is  an  extraordinarily  flexible  and 
powerful  modality  with  the  ability  to  report  on  multiple 
structural  and  functional  parameters  that  may  be  relevant  to 
lesion  response  assessment,  as  detailed  in  this  review. 

It  is  important  to  note  that  response  assessment  is  only 
one  of  two  goals  for  imaging  during  NAT,  with  the  other 
being  assessment  of  residual  disease  as  an  adjunct  to  surgi¬ 
cal  planning.11-14  While  it  is  tempting  to  think  of  these  two 
objectives  as  one  and  the  same,  assessment  for  residual 
disease  fundamentally  seeks  to  depict  anatomy,  whereas 
response  assessment  seeks  to  evaluate  changes  in  tumor 
biology  and  may  do  so  with  a  variety  of  techniques  reporting 
on  either  anatomical  or  functional  changes.  The  distinction 
is  especially  important  for  M  Rl,  where  certain  tradeoffs  in 
how  images  are  acquired  (eg,  between  high  temporal  and 
high  spatial  resolution  techniques)  may  theoretically  render 
one  imaging  protocol  more  useful  for  response  assessment 
and  another  protocol  more  useful  for  demonstrati  ng  ti  ny  foci 
of  residual  tumor.  This  review  focuses  on  M  Rl  techniques 
for  response  assessment,  but  we  address  practical  issues  of 
protocol  design  in  a  later  section  on  clinical  translation. 

Currently  available  quantitative  MRI 
techniques  for  response  assessment 

A  number  of  MRI  techniques  are  already  available  for 
potential  deployment  into  both  clinical  and  research  set¬ 
tings  for  the  assessment  of  breast  cancer  response  to  NAT. 
These  methods  include  anatomical  measurement  of  lesion 
size  (unidimensional,  multidimensional,  and  volumetric), 
DCE-M  Rl,  DW-M  R I ,  and  M  RS.  In  this  section  we  review  the 
basic  theory  underlying  these  techniques  (including  methods 
of  quantitative  analysis),  provide  examples  of  relevant  clini¬ 
cal  and  research  applications,  and  discuss  opportunities  for 
future  development. 

Lesion  size  measurement 

Theory 

By  virtue  of  its  exquisite  soft  tissue  contrast  and  very  high 
spatial  resolution,  MRI  is  an  extremely  powerful  technique 
for  demonstrating  the  morphology  of  breast  lesions.  M  Rl  has 
therefore  been  eval  uated  for  assessi  ng  breast  cancer  response 
to  therapy  using  changes  in  lesion  size  as  the  primary  mea¬ 
surement  variable.  M  R I -based  evaluation  of  lesion  size  relies 
al  most  exd  usi  vely  on  "contrast-enhanced"  i  magi  ng  ie,  i  mag- 
i  ng  fol  I  ow  i  ng  i  ntravenous  i  nj  ecti  on  of  a  paramagneti  c  contrast 


agent  (typically  a  gadolinium  chelate)  so  as  to  maximize  the 
conspicuity  of  an  enhancing  breast  lesion  against  a  back¬ 
ground  of  normal  breast  parenchyma.15  The  literature  has 
incorporated  uni  dimensional,  bi  dimensional,  and  volumetric 
(ie,  three-dimensional)  tumor  measurements  (Figure  1). 

Applications 

Evaluations  of  M  Rl-based  lesion  size  measurement  in  the 
setting  of  NAT  have  focused  on  two  clinical  scenarios: 

(1)  evaluation  for  residual  disease  at  the  end  of  NAT  and 

(2)  early  assessment  of  treatment  response.  With  regard  to 
evaluation  for  residual  disease,  a  significant  literature  has 
demonstrated  the  effectiveness  of  M  Rl-based  lesion  size 
measurement  in  predicting  results  at  surgical  pathology.5’16-22 
A  braham  et  al 19  showed  that  contrast-enhanced  M  Rl  at  the  end 
of  NAT  predicted  pathologically  determined  residual  diseasein 
97  percent  of  cases,  outperforming  both  physical  examination 
and  mammography.19  Balu-M  aestro  et  al 23  compared  physical 
examination,  mammography,  ultrasound,  and  M  Rl  and  found 
that  M  Rl  was  the  most  reliable  for  evaluating  residual  tumor 
size  after  NAT.  Londero  et  al22  found  that  M  Rl  after  NAT 
identified  residual  disease  better  than  mammography  and 
also  showed  that  M  Rl-measured  tumor  diameters  correlated 
better  with  pathologic  diameters  than  either  mammography  or 
sonography.  Chou  et  al 16  found  that  tumor  volume  after  NAT 
(calculated  using  a  segmentation  algorithm)  was  highly  cor¬ 
related  with  histopathologic  estimation  of  tissue  volume.16 

With  regard  to  early  assessment  of  treatment  response, 
a  number  of  studies  have  evaluated  the  ability  of  early 
changes  in  M  Rl-based  lesion  size  measurement  to  pre¬ 
dict  eventual  response.24-26  Martincich  et  al25  found  that  a 
minimum  65  percent  reduction  in  tumor  volume  after  two 
cycles  of  NAT  was  associated  with  an  eventual  major  his- 
topathological  response.  Padhani  et  al27  demonstrated  that 
a  decrease  in  bidimensional  tumor  area  after  two  cycles  of 
NAT  predicted  pathological  response.  Loo  et  al28  showed 
that  a  change  in  the  longest  diameter  of  enhancing  tumor 


Figure  I  Lesion  size  measurement  by  MRI.  (A)  Unidimensional  measurement  of 
tumor  long  axis  diameter.  (B)  Bidimensional  measurement  of  tumor  long  and  short 
axis  diameters.  (C)  Three-dimensional  measurement  of  tumor  volume. 
Abbreviation:  MRI,  magnetic  resonance  imaging. 
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predicted  pathologic  response  after  two  cycles  of  NAT.28 
The  recently  published  I  SPY-1  trial  showed  that  volumetric 
tumor  measurement  at  mid-treatment  outperformed  clinical 
assessment  in  predicting  pathologic  response.29 

Future  directions 

Although  M  Rl  techniques  for  lesion  size  measurement  are 
mature  and  have  been  shown  to  predict  changes  in  true  lesion 
size,  MRI-based  size  measurement  has  not  been  widely 
incorporated  into  the  clinical  or  research  settings  for  NAT 
response  assessment.  In  the  clinical  setting,  this  is  probably 
due  to  the  prevailing  opinion  that  physical  examination  pro¬ 
vides  adequate  response  assessment  in  the  patient  population 
for  whom  NAT  is  currently  indicated,  ie,  patients  with  large 
tumors.  In  the  research  setting,  RECIST-based  response 
biomarkers  may  be  incorporated  as  secondary  endpoints, 
but  current  trials  invariably  rely  on  pathologic  response  as 
the  primary  endpoint  for  determining  antitumor  efficacy. 
The  latter  remains  true  despite  some  preliminary  studies 
correlating  MRI-based  size  measurement  changes  with 
long-term  outcomes.30 

In  the  future,  with  the  anticipated  increased  use  of  NAT 
in  patients  with  smaller  tumors  that  are  poorly  assessed 
by  palpation,  MRI-based  lesion  size  measurement  may 
become  more  integrated  into  clinical  algorithms,  but  there 
remain  several  potential  problems  with  lesion  size-based 
response  assessment.  First,  many  breast  cancers  tend  to 
have  infiltrative,  irregular,  and/or  multifocal  growth,  creat¬ 
ing  problems  for  reproducible  size  measurement.  Second, 
as  mentioned  previously,  some  new  and  emerging  treat¬ 
ment  agents  have  cytostatic  rather  than  cytotoxic  effect 
and  may  produce  delayed  or  attenuated  lesion  shrinkage 
such  that  a  size-based  response  assessment  methodology 
may  not  fully  capture  relevant  changes  in  tumor  biology. 
Third,  when  treated,  breast  cancers  can  exhibit  a  variety 
of  different  morphological  changes  including  shrinkage 
from  the  outside,  melting  from  the  inside,  and  irregular 
internal  ("Swiss  cheese")  liquefaction;  these  different  pat¬ 
terns  would  be  poorly  assessed  by  methodologies  focusing 
exclusively  on  unidimensional,  bidimensional,  or  even 
volumetric  measurement. 

The  biggest  long-term  challenge  for  lesion  size 
measurement-based  techniques  will  therefore  probably 
involve  integration  with  other  techniques  focused  more  on 
assessing  functional  changes  in  tumor  biology,  as  described 
in  subsequent  sections.  Integration  of  different  techniques 
into  a  comprehensive,  multi-purpose  M  Rl  examination  may 
be  challenging,  and  will  also  be  described  later. 


Dynamic  contrast-enhanced 
MRI  (DCE-MRI) 

Theory 

DCE-MRI  involves  the  rapid,  sequential  acquisition  of 
Tj-weighted  images  before  and  after  the  injection  of  a  paramag¬ 
netic  contrast  agent.31  A  s  the  contrast  agent  perfuses  or  diffuses 
into  a  voxel  or  region  of  interest  (ROI),  it  shortens  the  native 
magnetic  relaxation  times  of  the  tissue  as  determined  by  the 
local  concentration  of  contrast.  W  hen  the  contrast  agent  leaves 
the  voxel,  the  relaxation  times  increase  toward  their  baseline 
value  at  a  rate  determined  by  local  tissue  characteristics.  Each 
voxel  thus  yields  a  signal  intensity  timecourse  that  can  be  ana¬ 
lyzed  to  yield  estimates  of  tissue  vascularity  parameters  includ¬ 
ing  perfusion,  permeability,  and  tissue  volume  fractions. 

Both  semiquantitative  and  fully  quantitative  methods 
have  been  developed  to  perform  DCE-MRI  analyses. 
Semiquantitative  methods  include  calculations  of  the  sig¬ 
nal  enhancement  ratio  (SER)  and  the  initial  area  under  the 
enhancement  curve  (iAUC). 32-36  The  SER  method,  used  in 
the  ISPY-1  trial,  employs  the  following  calculation: 

(1) 

“*2  -50 

where  S0  represents  the  signal  intensity  within  the  lesion 
before  contrast  administration,  S:  represents  the  signal  inten¬ 
sity  early  after  contrast  injection,  and  S2  represents  the  late 
postcontrast  signal  intensity.32  Given  a  map  of  SER  values, 
various  derivative  parameters  are  accessible  including  the 
SER  total  tumor  volume,  SER  partial  tumor  volume,  SER 
washout  tumor  vol  ume,  and  the  peak  SER. 

Fully  quantitative  DCE-M  Rl  analysis  typically  involves 
the  application  of  various  pharmacokinetic  equations  to 
model  the  movement  of  contrast  agent  molecules  between 
tissue  and  blood  vessels  over  time.37  Two  main  approaches 
have  been  deployed  in  breast  cancer:  the  standard Tofts-Kety 
(TK)  and  the  extended  Tofts-Kety  (ETK)  models.38  The 
standardTK  model  is  summarized  by: 

Ct(T)  =  K*ans  ■  jCp (t)  •  e-(KtarS/v^ '  (T_t)dt,  (2) 
o 

whereKtransisthe  volume  transfer  constant,  v  istheextravas- 
cular  extracellular  vol  ume  fracti  on,  Ct(T)  is  the  concentration 
of  contrast  in  the  tissue  of  interest,  and  Cp(t )  is  the  concentra¬ 
tion  of  contrast  in  blood  plasma  (also  known  as  the  arterial 
input  function).  The  ETK  model  incorporates  the  blood 
plasma  volume  fraction,  vp,  as  follows: 
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Ct(T)  =  /C^JQTCp(t)©<p(-(^ira'7\4)  (T-t))dt 
+  vPcP(t)- 

Both  theTK  and  ETK  models  assume  that  all  water 
compartments  within  tissue  are  well  mixed  (ie,  areatthefast 
exchange  I  i  mi t  of  the  nuclear  magneti c  resonance  ti  me  seal e) 
so  that  M  R I  signal  change  is  completely  described  by  asingle 
relaxation  rate  constant. 38  After  measurement  (or  estimation) 
of  Ct(t)  and  C  (t)  on  a  voxel  or  ROI  level,  Equations  2  and/or  3 
are  used  to  return  estimates  of  Ktrans,  ve,  and  v  ,  along  with  the 
derived  efflux  rate  constant  kep  (kep  =  Klrans/vJ.  The  working 
hypothesis  for  fully  quantitative  DCE-M  Rl  in  the  setting 
of  NAT  is  that  observed  changes  in  these  pharmacokinetic 
parameters  can  predict  treatment  response,  as  illustrated 
in  Figure  2.  For  a  more  extensive  discussion  of  DCE-M  Rl 
methods  in  oncology,  the  interested  reader  is  referred  to  the 
review  by  Yankeelov  and  Gore.31 

Applications 

Both  semiquantitativeand  fully  quantitative  DCE-M  Rl  tech¬ 
niques  have  been  evaluated  for  assessing  treatment  response 
in  breast  cancer  patients  undergoing  NAT.  The  best  known 
example  of  semi  quantitative  DCE-M  Rl  in  this  setting  is 
the  l-SPY  1  trial,29  in  which  patients  underwent  contrast- 
enhanced  M  Rl  at  multiple  time  points  before,  during,  and 
after  NAT;  the  authors  found  that  midtreatment  change  in 
SER  predicted  pathologic  response  with  an  area  under  the 
curve  (AUC)  of  0.71,  higher  than  clinical  assessment  (0.63) 
but  si  i  ghtly  I  ower  than  change  i  n  tumor  vol  ume  measurement 
(0.72).29  Other  investigations  of  semiquantitative  DCE-M  Rl 
include  employment  of  iAUC  for  response  assessment.  For 
example,  Tateishi  et  al36  used  the  percentage  area  under  the 
time-intensity  curve  (%AUC)  to  predict  treatment  response 
and  reported  that  although  the  sensitivity  of  %AUC  (50%) 


Figure  2  Fully  quantitative  DCE-MRI  analysis  in  a  breast  cancer  patient  undergoing 
NAT.  (A)  Pretreatment  DCE-MRI  analysis  yields  a  baseline  calculated  mean  tumor  Ktrans 
value  of  0.3  min-1.  (B)  DCE-MRI  analysis  after  one  cycle  of  NAT  yields  a  calculated  mean 
tumor  Ktrans  value  of  0.2  min-1.  (C)  Imaging  after  completion  of  NAT  shows  that  the 
lesion  is  no  longer  visible;  at  surgery,  the  patient  had  a  pathologic  complete  response. 
Note:  Ongoing  studies  are  investigating  whether  early  changes  in  mean  tumor  Ktrans 
can  reliably  differentiate  pathologic  responders  from  nonresponders. 
Abbreviations:  DCE-MRI,  dynamic  contrast-enhanced  magnetic  resonance 
imaging;  NAT,  neoadjuvant  therapy;  Ktrans,  volume  transfer  constant. 


was  not  acceptable,  the  specificity  of  %AUC  of  95.2%  was 
sufficiently  high  to  predict  pathologic  complete  response.36 
Our  group  has  shown  that  semiquantitative  analysis  of  high 
temporal  resolution  DCE-MRI  data  can  predict  pathologic 
response  after  a  single  cycle  of  NAT  (unpublished  data). 

Quantitative  analysis  of  DCE-M  Rl  data  has  been  shownto 
assess39  and  predict  treatment  response27  40  and  has  also  been 
shown  to  correlate  with  5-year  survival.41-44  Padhani  et  al27 
found  that  after  two  cycles  of  NAT,  both  change  in  tumor 
size  and  change  in  K trans  range  on  histogram  analysis  were 
equally  able  to  predict  pathologic  response.  Ah-See  et  al40 
analyzed  multiple  quantitative  DCE-M  Rl  parameters  and 
reported  that  change  i  n  K trans  was  the  best  predi  ctor  of  patho- 
logic  nonresponse.40  O  ur  group  found  that  not  only  the  mean 
but  also  the  standard  deviation  of  Ktans  as  estimated  by  the 
TK  and  ETK  models,  as  well  as  v  ,  can  separate  complete 
pathologic  responders  from  nonresponders  after  a  single 
cycle  of  NAT.45 

Diffusion-weighted  MRI  (DW-MRI) 

Theory 

Like  conventional  MRI,  DW-MRI  records  signals  from 
mobile  water  molecules  within  tissues,  but  in  DW-M  Rl  the 
contrast  reflects  the  distance  water  molecules  can  migrate 
or  "diffuse"  from  their  original  spatial  position  over  a  short 
time  interval  due  to  random,  thermally-induced  motion  (ie, 
Brownian  motion).  DW-M  Rl  exploits  applied  gradients  of  the 
main  magnetic  field  that  allow  for  localization  and  calculation 
of  the  microscopic  diffusion  of  water  molecules.  By  acquir¬ 
ing  two  or  more  images  with  different  degrees  of  "diffusion 
weighting"  (obtained  by  applying  the  diffusion-sensitizing 
gradients  with  different  amplitudes  on  successive  image 
acquisitions),  an  estimate  of  the  amount  of  molecular  water 
diffusion,  termed  the  apparent  diffusion  coefficient  (ADC), 
can  be  calculated  at  each  voxel  from 

S  =  S„exp(-b-ADC),  (4) 

where  5  is  the  signal  intensity  recorded  with  application  of 
a  diffusion-sensitizing  gradient,  S„  is  the  signal  intensity 
with  no  diffusion-sensitizing  gradient,  and  b  is  a  composite 
variable  reflecting  various  acquisition  parameters  (including 
the  strength  of  the  gradient  pulse,  duration  of  the  pulse,  and 
interval  between  pulses).46  ADC  values  from  successive 
voxels  can  then  be  aggregated  to  produce  a  map  of  ADC 
values  over  the  volume  of  interest  (Figure  3). 

In  general,  experimentally-measured  ADC  values  are 
lower  in  organized  tissues  than  in  free  solution  because 
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Figure  3  DW-MRI  in  a  breast  cancer  patient  undergoing  NAT.  (A)  On  a  pretreatment 
image  with  no  diffusion  gradient  (ie,  b  =  0  s/mm2),  the  tumor  is  difficult  to  distinguish 
from  background  normal  parenchyma.  (B)  Pretreatment  diffusion-weighted  image 
(b  =  660  s/mm2)  demonstrates  subtle  patchy  increased  signal  in  the  deep  upper 
breast,  corresponding  to  an  infiltrative  tumor.  (C)  Pretreatment  quantitative  ADC 
map,  with  color-coded  voxels  corresponding  to  tissue  ADC.  The  tumor  region  is 
outlined  in  white.  (D)  ADC  map  derived  from  DW-MRI  after  one  cycle  of  NAT;  the 
tumor  volume  (again  outlined  in  white)  has  markedly  decreased. 

Note:  This  patient  went  on  to  have  a  complete  pathologic  response. 
Abbreviations:  DW-MRI,  diffusion-weighted  magnetic  resonance  imaging;  NAT, 
neoadjuvant  therapy;  ADC,  apparent  diffusion  coefficient. 

various  structures  including  cell  membranes  and  intracel¬ 
lular  organelles  tend  to  restrict  or  hinder  the  free  movement 
of  water.47'48  Moreover,  cancerous  tissues  often  show  sig¬ 
nificantly  reduced  A  DC  values  when  compared  with  healthy 
tissues,  a  finding  typically  attributed  to  the  increased  cell 
density  of  many  malignancies.49  With  treatment,  intratumoral 
ADC  values  often  rise,  presumably  because  of  decreases  in 
cell  density  consequent  to  apoptosis  and  cell  death,  with 
concomitant  disruption  of  cell  membranes  allowing  water 
molecules  to  diffuse  more  freely.  This  basic  paradigm  -  low 
tumor  A  DC  values  before  treatment,  followed  by  rising  tumor 
ADC  values  with  treatment  -  provides  the  basic  model  for 
DW-M  R I  as  a  response  assessment  technique.  Importantly, 
it  has  been  shown  that  rising  ADC  values  can  occur  quite 
early  during  treatment,50’52  thus  providing  the  motivation  for 
studying  DW-M  Rl  as  an  early  response  biomarker.  For  a  more 
extensive  discussion,  the  reader  is  referred  to  the  review  by 
A  rl  i  nghaus  and  Y ankeel ov.53 

Applications 

Several  studies  have  correlated  changes  in  ADC  on  DW- 
MRI  with  treatment  response  in  breast  cancer  patients 


undergoing  NAT.  Early  investigations  demonstrated  increases 
in  meantumorADC  foil  owing  chemotherapy39  and  correlated 
mean  tumor  A  DC  increases  with  radiological  response,54'55 
while  subsequent  studies  examined  the  relationship  between 
mean  tumor  ADC  changes  and  pathological  response.  In 
general,  mean  tumor  ADC  has  been  found  to  increase  after 
NAT  in  both  pathologic  responders  and  nonresponders  but  to 
increase  more  for  pathologic  responders.56-59  Sharma  et  al52 
showed  that  at  the  end  of  therapy,  DW-M  R I  had  a  lower  sensi¬ 
tivity  but  a  higher  specificity  than  lesion  size  measurement  for 
differentiating  pathologic  responders  from  nonresponders. 

M  ore  recent  studies  have  looked  into  changes  in  ADC 
early  in  the  treatment  course  and  have  begun  to  investigate 
whether  DW-M  Rl  performed  early  in  treatment  can  be  used 
to  predict  response.  Pickles  etal51  demonstrated  a  significant 
increase  in  mean  tumor  ADC  after  a  single  cycle  of  NAT 
and  showed  that  mean  tumor  A  DC  began  to  change  earlier 
than  tumor  longest  diameter.  I  n  a  study  of  15  patients,  mean 
tumor  A  DC  increased  significantly  ( P  =0.008)  by  11%  after 
a  single  cycle  of  NAT.60  Li  etal50  reported  that  after  one  cycle 
of  NAT,  mean  tumorADC  significantly  increased  by  24%  in 
patients  who  went  on  to  have  a  complete  or  partial  response 
and  did  notchangein  patients  who  ultimately  exhibited  stable 
or  progressive  disease. 

Proton  magnetic  resonance 
spectroscopy  (MRS) 

Theory 

In  contrast  to  conventional  M  Rl,  which  generates  anatomi¬ 
cal  images  using  signal  primarily  arising  from  water,  M  RS 
provides  information  on  the  concentrations  of  different 
metabolites  in  tissue.  This  technique  can  be  exploited  to 
detect  the  altered  metabolic  signatures  of  cancer  cells.  For 
example,  many  malignancies  demonstrate  elevated  levels 
of  choline  and  lactate,  the  former  due  to  high  rates  of  mem¬ 
brane  turnover  and  the  latter  due  to  utilization  of  anaerobic 
glycolysis.  Choline  is  present  in  less  than  one  mi  1 1  i  molar 
concentrations  in  normal  breast  tissue  but  is  significantly 
el  evated  i  n  mal  i  gnant  breast  tumors6162  due  to  chol  i  ne  ki  nase 
overexpression  driven  by  H  IF-la63'64 

MRS  is  more  challenging  in  the  breast  than  in  other 
organs  due  to  the  large  amount  of  signal  from  lipid  and  to 
increased  magnetic  susceptibility  at  air-tissue  interfaces, 
although  the  latter  is  less  relevant  in  single-voxel  M  RS  of 
the  breast.  Despite  the  lower  spectroscopic  resolution  and 
signal-to-noise  ratio  of  breast  M  RS  compared  with  M  RS  in 
other  organs,  measurements  of  choline  levels  and  water-fat 
(W-F)  ratios  are  feasible  in  the  breast  and  have  been  used  for 
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breast  cancer  diagnosis  as  well  as  for  monitoring  response 
to  treatment  during  NAT.65'72 

Applications 

A  recent  clinical  trial  reported  by  Kumar  etal67  demonstrated 
that  malignant  breast  tissues  have  elevated  W-F  levels  com¬ 
pared  with  controls,  and  that  breast  cancers  decreasing  in  size 
with  NAT  also  exhibited  decreasing  W-F  ratios.67  In  a  sepa¬ 
rate  cl  i  ni  cal  tri  al  reported  byTozaki  etal71  using  choline  MRS 
in  patients  undergoing  NAT,  the  reduction  rates  of  choline 
were  statistically  significantly  different  between  pathological 
responders  and  nonresponders  after  two  treatment  cycles, 
with  positive  and  negative  predictive  values  of  choline  MRS 
of  89%  and  100%,  respectively.  M  oreover,  it  was  found  that 
the  predictive  power  of  choline  MRS  was  greater  than  that 
of  volumetric  tumor  measurements.  Danishad  et  al73  found 
that  choline  si gnal-to-noise  ratio  may  be  useful  in  predicting 
tumor  response  to  NAT. 

Future  directions 

Although  relatively  high  predictive  power  has  been  demon¬ 
strated  i  n  prel  i  mi  nary  tri  al  s,  the  w  i  despread  cl  i  ni  cal  i  mpl  emen- 
tation  of  W-F  and  choline  M  RS  for  breast  cancer  treatment 
response  assessment  will  confront  several  challenges.  First, 
customized  data  postprocessing  and  internal/external  signal 
referencing  is  required  to  convert  raw  spectroscopic  data  into 
quantitative  information,  especially  for  clinical  analysis  of 
multi-voxel  M  RS  data.  Second,  because  the  majority  of  tissue 
metabol  i  tes  i  n  vi  vo  besi  des  I  i  pi  ds  are  present  at  millimolar  con¬ 
centrations,  the  M  RS  sampling  voxel  must  typically  be  very 
large  (1-8  cm3)  in  order  to  achieve  sufficient  signal-to-noise 
ratio  for  very  dilute  metabolites,  and  this  low  spatial  resolu¬ 
tion  translates  into  limited  ability  to  interrogate  small  tumors 
and  to  report  on  intralesional  heterogeneity.  Improvements 
on  the  low  spatial  resolution  of  conventional  M  RS  may  result 
from  emerging  methods  including  multivoxel  chemical  shift 
imaging61'70'74  and  hyperpolarized  M  R,  discussed  below. 

Emerging  quantitative  MRI 
techniques  for  response  assessment 

Several  additional  M  Rl  techniques  are  on  the  horizon  for 
possible  future  use  in  assessing  breast  cancer  response  to 
NAT.  Some  of  these  techniques  have  already  been  deployed 
in  other  disease  sites.  This  section  discusses  possible  appli¬ 
cations  of  MT,  CEST,  and  MR  elastography  as  treatment 
response  biomarkers.  For  each  technique,  we  briefly  discuss 
the  underlyi  ng  theory  as  wel  I  as  opportuni  ti  es  and  prel  i  mi  nary 
applications  in  breast  cancer. 


Magnetization  transfer  (MT)  MRI 

Theory 

M  T  i  s  a  method  for  detecti  ng  and  quantifyi  ng  the  protons  asso- 
ciated  with  tissue  macromolecules.  These  macromolecular 
protons  cannot  be  routinely  detected  with  conventional  M  Rl; 
MT  takes  advantage  of  the  communication  between  macro- 
molecular  protons  and  water  to  elicitan  indirect  measurement 
of  their  properties.  MT  is  typically  performed  by  applying  a 
preparatory  off-resonance  radiofrequency  pulse  to  selectively 
saturate  protons  associated  with  macromolecules  and  then 
measuring  the  attenuation  of  the  water  signal  that  occurs  as 
a  result  of  the  exchange  of  spin  information  (magnetization 
transfer)  between  the  saturated  macromolecular  protons  and 
free  water.  The  magnetization  transfer  itself  takes  place  via 
dipolar  coupling  and/or  chemical  exchange.75 

The  typical  method  for  characterizing  the  MT  effect  is  to 
compare  the  signal  intensity  between  a  nonsaturated  image 
(or  so-called  reference  acquisition,  designated  as  M  Toff)  and 
an  image  acquired  after  application  of  the  off-resonance 
radiofrequency  (RF)  pulse  ( MTon ).  The  magnitude  of  the 
saturation  is  proportional  to  the  quantity  of  saturated  mac¬ 
romolecular  protons  and  the  efficiency  of  exchange  (or 
exchange  rate)  with  free  water  protons  (Figure  4).  The  MT 
effect  is  thus  characterized  by  the  magnetization  transfer 
ratio  (MTR): 

MTR  =1-M  TJM  Toff,  (5) 

where  MTon  represents  the  signal  intensity  with  the  satura¬ 
tion  pulse  and  MToff  represents  the  signal  intensity  from  the 
reference  image.  The  M  TR  has  been  shown  to  be  related  to 
the  amount  of  macromolecular  protons  in  tissue.76  It  should 
be  noted  that  since  the  MTR  measurement  is  affected  by  the 
relaxation  and  exchange  rates,  the  field  strength,  and  the 
RF  irradiation  power,  it  is  only  a  semiquantitative  metric. 
A  quantitative  MT  approach  is  necessary  to  separate  the 


Figure  4  MT  results  from  a  healthy  volunteer.  (A)  MToff.  (B)  MTon.  (C)  MTR  map 
demonstrating  an  average  40%  reduction  in  signal  (ie,  MTR  =  0.4)  in  the  fibroglandular 
tissue  with  good  fat  suppression. 

Abbreviations:  MT,  magnetization  transfer;  MTon,  signal  intensity  with  the 
saturation  pulse;  MToff,  signal  intensity  from  the  reference  image;  MTR,  magnetization 
transfer  ratio. 
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contributions  from  MT  and  relaxation  effects,7778  and  the 
development  of  cl  i  ni  cal  ly  feasi  bl  e  quantitative  M  T  techni  ques 
is  an  active  area  of  research.79  For  a  more  extensive  discus¬ 
sion  of  M  T  methods,  the  interested  reader  is  referred  to  the 
review  by  Gochberg  and  Lepage.80 


Opportunities  and  preliminary  applications 
in  breast  cancer 

While  the  primary  application  of  MT  has  been  in  studies 
of  demyelinating  diseases,81-83  the  demonstration  of  MT's 
sensitivity  to  collagen  content  in  meningiomas84  has  led 
to  interest  in  MT  for  studying  changes  in  the  extracellular 
matrix  (ECM )  of  the  breast. The  EC  M  is  a  major  component 
of  the  fibroglandular  tissue  of  the  breast  and  is  comprised 
of  a  network  of  macromolecules,  including  collagen, 
fibronectin,  and  laminin.  The  ECM  is  known  to  play  a  role 
in  tumor  development  and  progression,85-87  and  the  role  of 
the  ECM  in  breast  cancer  has  gained  increasing  attention. 
For  example,  loachim  et  al88  reported  that  the  expression  of 
the  ECM  macromolecules  fibronectin,  collagen  type  IV,  and 
laminin  is  altered  in  breast  cancer;  Levantal  et  al89  recently 
demonstrated  that  cross-1  i nking  of  collagen  type  I  isi nvol ved 
in  the  stiffening  of  the  ECM  and  tumor  migration  in  breast 
cancer;  and  Kakkad  et  al90  reported  that  the  concentration  of 
collagen  type  I  is  reduced  in  the  hypoxic  tumor  environment. 
An  initial  application  of  MT  in  breast  cancer  showed  a 
significant  reduction  of  MTR  in  malignancies  compared 
with  benign  tumors,  thought  to  represent  a  reduction  of  the 
macromolecular  pool  due  to  increased  proteolytic  activity.91 
In  the  NAT  setting,  it  is  hypothesized  that  changes  in  mac¬ 
romolecular  content  in  response  to  successful  treatment 
would  result  in  measurable  changes  in  MTR  values.  Also, 
as  tumors  respond  to  therapy  and  shrink,  the  tumor  cells  are 
typically  replaced  with  fibrotic  tissue,92-94  which  potentially 
could  lead  to  an  (ultimate)  increase  in  MTR  in  the  case  of 
successful  therapy. 

Chemical  exchange  saturation 
transfer  (CEST)  MRI 

Theory 

CEST  is  similar  to  MT,  but  rather  than  focusing  on 
macromolecules,  CEST  seeks  to  specifically  irradiate  tissue 
metabolites  such  as  amides,  amines,  and  hydroxyl  groups  that 
are  also  in  exchange  with  free  water.95-97  The  exchangeable 
protons  on  these  metabolites  have  chemical  shifts  that  are 
significantly  smaller  than  the  broad  macromolecular  pool 
and  exchange  at  significantly  slower  rates,  which  allows  for 


spectrally  selective  irradiation  via  application  of  an  off- res¬ 
onance  RF  pulse.98'99  Under  experimental  conditions,  CEST 
has  been  shown  to  discriminate  individual  tissue  metabolites 
with  high  specificity.  In  addition,  because  proton  chemical 
exchange  rates  are  pH  -dependent,  CEST  can  be  used  to  inter¬ 
rogate  for  changes  in  tissue  pH.100  Like  MT,  CEST  can  be 
performed  without  the  use  of  exogenous  contrast  agents. 

CEST  is  performed  through  application  of  a  spectrally- 
selective  saturation  pulse  prior  to  an  imaging  sequence,  as 
shown  in  Figure  5.  This  series  is  repeated  while  the  reso¬ 
nance  offset  of  the  saturation  pulse  is  swept  through  a  range 
of  frequencies,  typically  <  ±10  ppm.  The  saturation  will 
affect  specific  protons,  and  the  observed  water  signal  will  be 
attenuated  viadirect  chemical  exchange.100  CEST  results  are 
often  exami  ned  vi  a  a  z-spectrum,  a  pi  ot  of  the  si  gnal  i  ntensi  ty 
of  water  as  a  function  of  saturation  offset,  normalized  by 
the  signal  intensity  of  water  in  the  absence  of  saturation.101 
Example  z-spectra  are  shown  in  Figure  6,  which  depicts  the 
results  of  CEST  analysis  at  3 T  on  a  malignant  breast  tumor 
(black  line)  compared  with  healthy  fibroglandular  tissue 
(gray  line).  One  way  the  CEST  effect  can  be  characterized 
is  by  examining  the  amount  of  asymmetry  observed  in  the 
CEST  spectra,  with  the  effect  at  particular  offset  frequencies 
related  to  the  exchanging  protons  of  interest. 

Perhaps  the  most  widely  reported  CEST  effect  is  derived 
from  the  exchange  of  amide  protons  on  the  backbone  of 
proteins  and  peptides.  In  this  case,  the  CEST  effect  of  the 
amide  proton  resonancecan  be  characterized  by  the  so-called 
proton  transfer  ratio  (PTR),  a  measure  of  the  asymmetry 
of  the  z-spectrum  about  the  water  frequency  ( CESTasym )  at 
A co- 3.5  ppm: 

_  S{-Aa>) -S(Aai) 

LC-)'asym  <-  '  '°' 

-’o 

where  S(-Aa>)  and  S(Aa>)  are  the  signal  intensities  with  the 
saturation  at±A<w,  andS0  is  the  signal  intensity  in  the  absence 
of  saturation.  This  calculation  negates  the  confounding  effects 
of  direct  water  saturation,  which  are  symmetric  about  the 
center  of  the  z-spectrum. 102  The  measured  PTR  can  be  affected 
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> 
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acquisition 


Figure  5  General  pulse  sequence  diagram  for  a  CEST  MRI  experiment. 

Note:  RF  irradiation  for  a  time  ts  with  an  amplitude  of  precedes  the  excitation  (a 
degrees)  and  image  acquisition. 

Abbreviations:  CEST  MRI,  chemical  exchange  saturation  transfer  magnetic 
resonance  imaging;  RF,  radiofrequency. 
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Figure  6  Example  z-spectra  arising  from  a  CEST  MRI  experiment  at  3  T. 

Note:  The  normalized  signal  (S/SQ)  is  shown  as  a  function  of  saturation  offset 
frequency  for  regions  of  interest  in  malignant  tumor  (black  line)  and  healthy 
fibroglandular  tissue  (gray  line). 

Abbreviations:  CEST  MRI,  chemical  exchange  saturation  transfer  magnetic 
resonance  imaging;  S,  signal  intensity  with  saturation;  S0,  signal  intensity  in  the 
absence  of  saturation. 


by  the  concentrati  on  of  exchangi  ng  protons  as  wel  I  as  the  pH , 
which  infl uences  the  proton  exchange  rate.  E xtensi ve  research 
is  underway  to  determine  the  underlying  mechanism  driving 
the  measured  CEST  effect  in  both  healthy  and  diseased  tis¬ 
sue.  Ongoing  research  is  also  investigating  innovative  pulse 
sequences103'104  and  alternative  quantification  strategies105-107 
that  may  minimizeasymmetric  magnetization  transfer  effects 
from  macromolecules,  which  can  confound  the  PTR  measure¬ 
ment.108  Fora  more  comprehensive  i  ntroduction  to  C  E  ST,  the 
reader  is  referred  to  the  review  by  Gochberg  and  Lepage.80 


Opportunities  and  preliminary 
applications  in  breast  cancer 
CEST  focused  at  the  amide  proton  resonance  (3.5  ppm) 
provides  information  on  the  amide  protons  of  protein/ 
peptide  backbones  and  has  been  used  to  study  tissues  where 
either  the  protein/peptide  concentration  or  the  pH  may  be 
altered.95'109-111  This  technique,  sometimes  termed  amide 
proton  transfer  (APT)  imaging,  is  thought  to  be  especially 
relevant  to  cancer  imaging  because  tumor  cells  may  accu¬ 
mulate  defective  protei  ns  at  a  hi  gher  rate  than  normal  and/or 
experience  alterations  in  pH  due  to  hypoxia.112  A  PT-CEST 
has  been  applied  to  brain,110  prostate,113  and  breast  tumors106 
and  has  been  used  to  differentiate  cellular  protein  content 
between  tumor  and  healthy  cells.114 

0  ur  group  is  actively  developi  ng  A  PT-CEST  i  maging  as  a 
potential  early  treatment  response  bi  omarker  for  breast  cancer 
NAT.106115  Figure  7  depicts  the  measured  APT  overlaid  on 
an  anatomical  image  for  three  breast  cancer  patients  before 
(left  column)  and  after  (right column)  onecycleof  NAT.The 


Figure  7  Amide  proton  transfer  (APT)  maps  derived  from  CEST  MRI  in  breast 
cancer  patients  undergoing  NAT.  Baseline  (pretreatment)  images  are  presented  on 
the  left,  and  images  after  one  cycle  of  NAT  are  presented  on  the  right.  (A):  patient 
with  complete  response  after  one  cycle  of  therapy  (27%  decrease  in  measured  APT 
from  baseline).  (B):  patient  with  partial  response  (49%  increase  in  measured  APT). 
(C):  patient  with  progressive  disease  (78%  increase  in  measured  APT). 
Abbreviations:  CEST  MRI,  chemical  exchange  saturation  transfer  magnetic 
resonance  imaging;  NAT,  neoadjuvant  therapy. 


top  row  shows  imaging  data  from  a  complete  responder  (27% 
decrease  in  measured  A  PT  from  baseline),  the  middle  row 
from  a  partial  responder  (49%  increase  in  measured  APT), 
and  the  bottom  row  from  a  patient  with  progressive  disease 
(78%  increase  in  measured  APT).  These  preliminary  results 
demonstrate  the  potential  sensitivity  of  A  PT-CEST  to  the 
molecular  changes  occurring  early  during  treatment. 

In  addition  to  A  PT-CEST,  there  is  potential  to  examine 
breast  cancer  by  deploying  CEST  to  evaluate  tissue  gly- 
cosaminoglycan  content  (ie,  gagCEST).  M  any  different  cell 
surface  and  matrix  proteoglycan  core  proteins  are  expressed 
i  n  the  mammary  gl ands.  T  he  level  of  expressi on  of  these  core 
protei  ns,  the  structure  of  thei  r  glycosami  noglycan  chai  ns,  and 
their  degradation  are  regulated  by  many  of  the  same  effectors 
that  control  development  and  function. 116  Loss  or  overexpres¬ 
sion  of  proteoglycans  in  carcinoma  cells  has  been  associated 
with  malignant  progression117'113  and  has  correlated  with  poor 
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prognosis,119  leading  to  the  hypothesis  that  gagC  EST  may  be 
developed  as  a  treatment  response  biomarker  in  the  future. 
T  he  appl  i  cati  on  of  C  E  ST  targeti  ng  si  gnatures  of  active  tumors 
bri  ngs  potential  for  noni  nvasi ve  mol ecul ar  i  magi  ng  that  could 
be  predictive  of  prognosis. 

Magnetic  resonance  elastography  (MRE) 

Theory 

MRE  is  based  on  use  of  the  elastic  properties  of  tissue  as 
an  imaging  contrast  mechanism.  The  general  concept  of 
elastography,  realized  first  in  ultrasound120  and  later  devel¬ 
oped  in  MRI,121  involves  the  use  of  imaging  to  measure 
tissue  response  to  applied  physical  deformation.  Sometimes 
described  asaform  of  "21stcentury  palpation,”  elastography 
allows  for  generation  of  tissue  elasticity  maps  providing 
spatial  visualization  and  quantification  of  the  distribution  of 
elastic  properties  within  an  object.  Generally,  elastography 
methods  employ  the  simplifying  assumption  of  a  linear 
elastic  isotropic  constitutive  model  where  mechanical  equi¬ 
librium  is  governed  by: 

V  •  GVO+V^^(V  •  D)  +  F  =-(/w-/W9. 

where  G  is  the  shear  modulus,  u  is  the  displacement  vec¬ 
tor,  v  is  Poisson's  ratio,  F  are  body  forces,  p  is  density,  and 
g  is  the  gravitational  constant  vector.  Given  a  measured 
tissue  displacement  field  along  with  appropriate  assump¬ 
tions,  the  above  equation  can  be  used  to  reconstruct  the 
spatial  distribution  of  shear  modulus  (orYoung's  modulus, 
E  =  2G[1+ v]). 

MRE  methods  can  be  broadly  classified  as  either  dynamic 
or  static.  With  dynamic  excitation  MRE,  shear  waves  are 
applied  to  an  area  of  interest  by  piezoelectric  or  pneumatic 
sources;  oscillating  motion-sensitized  gradient  sequences 
are  synchronized  to  the  externally-applied  excitation,  and 
induced  three-dimensional  tissue  motion  is  recorded  with 
phase-contrast  imaging.121-124  Static  MRE  methods125-128 
reconstruct  elastic  properties  of  tissue  by  incorporating 
a  biomechanical  finite  element  model  into  the  nonrigid 
registration  of  images  acquired  under  different  static  load¬ 
ing  environments;  the  application  of  this  technique  is  not 
unique  to  M  Rl  modality,  and  in  fact  this  methodology  is  also 
referred  to  as  modality-independent  elastography  (MIE)128 
due  to  its  foundation  in  image  processing.  A II  M  RE  methods 
require  specially  designed  hardware  to  couple  the  deforma¬ 
tion  source  (either  static  or  dynamic)  to  the  area  of  interest. 
Challenges  to  successful  application  of  MRE  include 


synchronization  of  the  dynamic  deformation  source  to  the 
motion-sensitized  M  R  gradient  and  often  complex  biome¬ 
chanical  model-driven  inversions  of  the  deformation  field  to 
reconstruct  tissue  elasticity  maps.  For  a  more  extensive  review 
of  M  R-based  elastography  techniques,  the  interested  reader  is 
referred  to  M  ariappan  et  al129  and  M  anduca  et  al.130 

Opportunities  and  preliminary  applications 
in  breast  cancer 

The  rationale  for  exploring  MRE  as  a  possible  breast  cancer 
treatment  response  biomarker  derives  from  observations  that 
cancer  progression  is  accompanied  by  extracellular  matrix 
remodeling  and  increasing  mechanical  stiffness.89  M  RE  has 
been  deployed  in  liver  tumor  assessment  and  has  had  promis¬ 
ing  preliminary  success  in  distinguishing  between  benign  and 
malignant  lesions,  with  malignant  tumors  exhibiting  a  sig- 
nificant  (>3-4  fold)  stiffness  increase  over  benign  tumors.131 
While  there  are  very  few  studies  of  M  R-based  elastography 
in  breast  cancer,  preliminary  results  show  initial  promise  in 
quantifying  in  vivo  stiffness  of  breast  tumors  (as  well  as  dif¬ 
ferent  ati on  of  fi brogl andul ar  and  adi  pose  ti ssue)132  and  have 
demonstrated  the  ability  to  improve  diagnostic  sensitivity  and 
specificity  over  more  traditional  breast  M  Rl.133 

Our  group  is  actively  developing  M  RE  as  a  mechanical 
property  biomarker  in  breast  cancer,  and  our  novel  MIE 
method  has  been  under  preliminary  investigation  in  breast 
cancer.126’127  Recent  breakthroughs  involving  translation  and 
automation127  have  generated  promising  preliminary  results 
towards  the  ability  to  spatially  reconstruct  tissue  elastic 
mechanical  properties  in  breast  cancer.  Figure  8  depicts  the 
results  of  M  I E  analysis  on  a  breast  cancer  patient. 


Figure  8  Static  MRE  (MIE)  results  from  a  breast  cancer  patient.  (A)  undeformed 
image  volume.  (B)  deformed  image  volume.  (C)  undeformed  central  slice. 
(D)  deformed  central  slice.  (E)  reconstructed  tissue  elasticity  map. 
Abbreviations:  MRE,  magnetic  resonance  elastography;  MIE,  modality-independent 
elastography. 
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M  onitoring  the  changes  in  mechanical  properties  along 
the  NAT  time  course  may  provide  a  means  of  assessing 
response  to  therapy.  Additionally,  quantitative  mechanical 
property  information  can  also  be  used  to  inform  patient- 
specific  reaction-diffusion  tumor  growth  models  that  are 
mechanically  coupled  to  constrained  tumor  cell  diffusion.134 
The  growth  and  tumor  cell  diffusion  parameters  within  the 
tumorgrowth  modelscan  befitbetween theinitial  and  inter¬ 
mediate  time  points  and  then  projected  forward  to  the  final 
time  point  in  order  to  help  predict  patient  response.135'136 

Hyperpolarized  MR 

Theory 

Hyperpolarized  M  R  seeks  to  overcome  the  limitations  of 
conventional  M  RS  by  exploiting  exogenous  contrast  agents 
that  have  been  "hyperpolarized,”  ie,  a  large  proportion  of 
their  nuclear  spins  have  been  aligned  with  the  magnetic  field 
or  polarized.  Whereas  conventional  M  R  imaging  depends  on 
spins  that  have  been  polarized  on  the  order  of  a  few  parts  per 
million,  hyperpolarization  allows  for  polarizing  nuclear  spins 
to  nearly  unity.  The  resulting  several-orders-of-magnitude 
i ncrease i n  signal  i ntensity  translates  i nto  i ncreased  sensitivity 
for  detecting  metabolic  markers  of  cancer  such  as  lactate  and 
bicarbonate  that  may  be  present  at  millimolar  and  potentially 
submillimolar  concentrations137'138  and  that  may  not  be  ame¬ 
nable  to  conventional  M  RS.137'139  The  technique  may  also  be 
exploited  to  provide  better  spatial  resolution  and  significantly 
faster  examination  times  compared  with  conventional  proton 
MRS.  Long-lived  nuclear  spin  sites  are  typically  selected 
for  the  preparation  of  hyperpolarized  contrast  agents  that 
are  suitable  for  injection  into  living  organisms.  Carbon-13 
(13C)  sites  without  directly-attached  protons  are  most  com¬ 
monly  used,  due  to  their  long  T2  of  -  20-40  seconds  and  their 
abundancewithin  many  metabol  ically-relevant  molecules.  At 
present,  clinical  translation  of  hyperpolarized  MR  technol¬ 
ogy  is  limited  by  demanding  instrumentation  and  software 
requirements,  including  multinudear  M  Rl  scanner  capability, 
highly  specialized  RF  pulse  sequences  (tailored  to  the  RF  coil, 
magnetic  field  strength,  and  metabolic  contrast  agent),  and 
multinudear  RF  coils. 

Opportunities  and  preliminary  applications 
in  breast  cancer 

While  in  vivo  hyperpolarized  M  R  technology  is  less  than 
10  years  old,  it  has  rapidly  progressed  from  proof-of-principle 
studies  in  miceto  the  first  clinical  trial  in  humans,  including 
analyses  of  choline,  pyruvate,  fumarate,  and  bicarbonate 
as  potential  biomarkers.137'140-143  There  have  been  as  yet  no 


specific  clinical  applications  to  breast  cancer,  but  preliminary 
biomarker  studies  in  mice  and  in  human  prostate  cancer  sug¬ 
gest  this  technique  may  be  applicable  for  breast  cancer  NAT 
response  assessment  in  the  future.143 

Challenges  and  opportunities 
for  clinical  translation 

Several  important  challenges  must  be  addressed  in  work¬ 
ing  toward  the  translation  and  adoption  of  quantitative  M  Rl 
techniques.  First,  investigators  will  have  to  choose  meaningful 
and  clinically  relevant  statistical  methods  for  validating  novel 
M  Rl  methods.  Prior  studies  of  M  Rl  as  a  response  assessment 
tool  for  NAT  in  breast  cancer  have  chosen  a  variety  of  clinical 
outcome  variables  -  including  palpation,  imaging-based  size 
measurement  change,  and  pathologic  response  -  against  which  to 
evaluate  imaging  as  a  response  biomarker.  Pathologic  response 
is  the  most  clinically  relevant  of  these  outcome  variables,  hav¬ 
ing  been  established  as  an  independent  prognostic  marker  for 
overall  survival  in  breast  cancer,3  yet  researchers  attempting  to 
validate  their  methods  against  pathologic  response  will  have 
to  contend  with  the  multiple  different  definitions  of  pathologic 
response  now  circulating  in  the  breast  cancer  literature  as  well 
as  the  notion  raised  in  certain  recent  studies  that  pathologic 
response  may  be  a  suitabl  e  surrogate  endpoi  nt  for  some  but  not 
all  breast  cancer  subtypes.144  Widespread  incorporation  of  quan¬ 
titative  M  Rl  response  biomarkers  in  the  clinical  and  research 
settings  will  probably  require  validation  not  just  as  predictors 
of  pathologic  response,  but  as  predictors  of  long-term  outcomes 
including  recurrence  rates  and  overall  survival. 

A  second  challenge  involves  ensuring  standardization  and 
reproducibility  of  M  Rl  methods.  M  R  imaging  is  a  complex 
undertaking,  with  images  and  quantitative  measurements 
dependent  upon  multi  pie  precise  software  setti  ngs  and  hardware 
configurations.  Results  of  a  controlled  experiment  on  a  single 
scanner  in  a  laboratory  environment  may  be  quite  difficult 
to  reproduce  across  multiple  imaging  sites.  Researchers  will 
have  to  work  with  vendors  to  increase  reproducibility  across 
platforms,  and  vendors  will  likely  be  called  upon  to  provide 
increased  transparency  with  regard  to  proprietary  hardware  and 
software  architectures.  A  handful  of  government- industry  con¬ 
sortia  including  the  Radiological  Society  of  North  America's 
Quantitative  Imaging  Biomarkers  Alliance  are  attempting  to 
f  aci  I  itate  these  devel  opments.  We  may  al  so  see  conti  nued  i  ncor- 
poration  of  quantitative  MRI  techniques  into  commercially 
available  computer-assisted  detection  software. 

A  third  challenge  for  clinical  translation  lies  in  design¬ 
ing  M  R  imaging  protocols  to  meet  the  various  objectives 
of  a  clinical  scan.  A  s  mentioned  earlier,  there  are  two  chief 
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goals  for  breast  imaging  during  NAT  -  response  assessment 
and  demonstration  of  residual  disease  -  and  these  goals 
may  call  for  different  M  Rl  approaches  that  may  be  mutu¬ 
ally  exclusive.  For  example,  basic  M  Rl  principles  dictate 
a  fundamental  tradeoff  between  spatial  and  temporal 
resolution-in  general,  high  spatial  resolution  images  require 
longer  acquisition  times  and  therefore  a  necessary  sacrifice 
of  temporal  resolution,  and  high  temporal  resolution  images, 
by  virtue  of  the  rapidity  with  which  they  are  acquired,  do 
not  provide  sufficient  time  to  gather  high  spatial  resolution 
data.  A  scan  tailored  for  demonstration  of  residual  disease 
may  call  for  high  spatial  resolution  imaging,  whilea  DCE- 
M  Rl  scan  tailored  for  early  response  assessment  may  call 
for  high  temporal  resolution  imaging  in  order  to  provide 
the  most  accurate  model  of  vascular  flow.  These  tradeoffs 
highlight  the  challenges  in  clinical  MRI  protocol  design, 
where  different  clinical  objectives  sometimes  compete  with 
each  other.  Potential  solutions  may  emerge  from  creating 
hybrid  protocols  incorporating  a  variety  of  imaging  tech¬ 
niques;  from  exploring  newer  methods,  such  as  compressed 
sensing,  that  may  make  the  above  tradeoffs  less  apparent; 
and/or  from  employing  different  scanning  protocols  for 
different  settings  such  as,  a  high  temporal  resolution  M  Rl 
protocol  early  in  the  course  of  therapy  and  a  high  spatial 
resolution  protocol  at  the  end  of  therapy  as  an  adjunct  to 
surgical  planning. 

Finally,  investigators  across  multiple  disciplines  will 
have  to  work  toward  overcoming  the  challenge  of  infor¬ 
mation  saturation  that  clinicians  will  inevitably  face  as 
the  methods  discussed  in  this  article  are  translated  into 
the  clinical  setting.  As  more  quantitative  M  Rl  metrics  are 
validated  as  useful  biomarkers,  we  may  see  development 
of  integrated  clinical  scoring  systems  that  synthesize 
multi  parametric  imaging  variables  with  other  quantitative 
clinical  parameters. 

Conclusion 

There  is  growing  demand  for  objective  and  standardized 
early-response  biomarkers  for  breast  cancer  NAT,  with  such 
biomarkers  to  be  used  as  an  adjunct  to  and  possibly  eventu¬ 
ally  a  replacement  for  pathologic  assessment  of  treatment 
response. The  imaging  research  community  is  responding  to 
this  demand  by  developing  and  testing  novel  approaches  in 
preclinical  models,  in  single  site  studies,  and  in  large  multi  site 
clinical  trials.  Quantitative  MRI  techniques,  including  the 
current  and  emerging  methods  reviewed  in  this  contribu¬ 
tion,  hold  great  potential  for  incorporation  into  both  clinical 
and  research  algorithms.  As  these  techniques  are  validated 


and  correlated  with  long-term  clinical  outcomes,  we  may 
witness  a  broad  transformation  in  the  use  of  imaging  with 
breast  cancer  NAT. 
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The  exquisite  NMR  spectral  sensitivity  and  negligible  reactivity 
of  hyperpolarized  xenon-129  (HP129Xe)  make  it  attractive  for  a 
number  of  magnetic  resonance  applications;  moreover,  HP129Xe 
embodies  an  alternative  to  rare  and  non-renewable  3  He.  However, 
the  ability  to  reliably  and  inexpensively  produce  large  quantities 
of  HP129Xe  with  sufficiently  high  129Xe  nuclear  spin  polarization 
(PXe)  remains  a  significant  challenge — particularly  at  high  Xe  den¬ 
sities.  We  present  results  from  our  open-source’  large-scale  (— 1 
L/hr)  129Xe  polarizer  for  clinical,  pre-clinical,  and  materials  NMR 
and  MRI  research.  Automated  and  comprised  mostly  of  off-the- 
shelf  components,  this  ‘hyperpolarizer’  is  designed  to  be  readily 
implementable  in  other  laboratories.  The  device  runs  with  high 
resonant  photon  flux  (up  to  200  W  at  the  Rb  D ^  line)  in  the  xenon- 
rich  regime  (up  to  1800  Torr  Xe  in  500  cc)  in  either  single-batch 
or  stopped-flow  mode,  negating  in  part  the  usual  requirement  of 
Xe  cryo-collection.  Excellent  agreement  is  observed  among  four 
independent  methods  utilized  to  measure  spin  polarization.  In¬ 
cell  PXe  values  of  —90%,  —57%,  —50%,  and  —30%  have  been 
measured  for  Xe  loadings  of  —300,  —500,  —760,  and  —1570  Torr, 
respectively.  Pxe  values  of  —41  %  and  —28%  (with  —760  and  — 1 545 
Torr  Xe  loadings)  have  been  measured  after  transfer  to  Tedlar 
bags  and  transport  to  a  clinical  3  T  scanner  for  MR  imaging  as  a 
prelude  to  in  vivo  experiments.  Long  in-bag’  129Xe  polarization 
decay  times  have  been  measured  (H- 38  min  and  —5.9  hr  at  —1.5 
mT  and  3  T,  respectively) — more  than  sufficient  for  a  variety  of 
applications. 


hyperpolarization  |  NMR&MRI  |  laser-polarized  xenon 


Owing  to  the  detection  sensitivity  provided  by  their  high,  non¬ 
equilibrium  nuclear  spin  polarization,  hyperpolarized  (HP)  noble 
gases  (e.g.  129Xe  and  3He)  are  utilized  in  a  number  of  NMR/MRI 
applications  (1,  2).  Human  lung  imaging  (e.g.,  Refs.  (3-7))  with 
HP129Xe  is  of  particular  interest.  Moreover,  xenon  is  soluble  in 
blood  (8),  other  tissues  (9),  and  many  biologically  compatible 
liquids  (10),  and  its  proclivity  for  interacting  with  other  substances 
and  its  wide  chemical  shift  range  make  HP129Xe  a  sensitive  NMR 
probe  of  molecular  and  material  surfaces  (1,  2,  11-13).  In  many 
NMR  and  MRI  applications  HP129Xe  can  replace  3He,  and  the 
relative  abundance  of  129Xe  can  greatly  reduce  the  impact  of  the 
world-wide  3  He  shortage  (14)  in  these  instances.  Despite  con¬ 
siderable  progress  (e.g.,  Refs.  (15-23)),  a  major  obstacle  towards 
implementing  HP129Xe  for  clinical  imaging  has  been  the  ability 
to  reliably  and  inexpensively  produce  large  quantities  of  HP129Xe 


with  high  polarization  (Pxe)[l].  HP  Xe  is  usually  created  via 
spin-exchange  optical  pumping  (SEOP  (24)),  a  process  whereby 
the  unpaired  electronic  spins  of  an  alkali  metal  vapor  (typically 
Rb)  are  polarized  via  depopulation  optical  pumping  with  circu¬ 
larly  polarized  laser1  light,  and  the  polarization  is  subsequently 
transferred  to  noble  gas  nuclear  spins  during  collisions.  It  is  gener¬ 
ally  anticipated  that  high  Px&  values  are  achievable  via  SEOP  only 
in  the  low  xenon-density  regime  (19,  25),  because:  ( a )  higher  Xe 
densities  increase  the  destruction  of  the  alkali  metal  polarization 
from  non-spin-conserving  collisions  at  a  rate  that  is  orders  of 
magnitude  worse  than  those  of  other  gases  like  N2  and  He  (26- 
28);  and  ( b )  higher  total  pressures  tend  to  quench  the  three-body 
van  der  Waals  contribution  to  Rb-Xe  spin  exchange — leaving  the 
less-efficient  two-body  term  to  dominate  (19,  24).  Most  large- 
scale  polarizers,  in  particular  all  that  are  available  commercially, 
operate  in  this  low-Xe  density  regime.  Applied  research  with  HP 
129Xe  is  severely  hampered  owing  to  a  lack  of  access  to  expensive 
proprietary  commercial  hyperpolarizers,  and  in  order  to  stimulate 
development  of  HP129Xe  applications,  we  present  here  a  low  cost, 
open-source  design. 

In  our  recent  work  (23,  29,  30)  exploring  batch-mode  and 
stopped-flow  (16)  Rb/Xe  SEOP  under  conditions  of  high  reso¬ 
nant  laser  flux,  an  unexpected,  inverse  relationship  was  found 
between  the  optimal  temperature  for  performing  Rb/Xe  SEOP 
and  the  in-cell  Xe  density.  This  effect  was  recently  exploited  to 
achieve  surprisingly  high  Pxe  values  at  high  in-cell  Xe  densities 
(e.g.,  -52%,  -31%,  -22%,  and  -11%  at  50,  300,  500,  and  2000 
Torr  Xe  in  a  75  cc  cell)  with  <  30  W  of  laser  power  (23).  Expanding 
these  small-scale  experiments  by  over  an  order  of  magnitude 
and  improving  PXe  is  a  scientific  and  engineering  leap  necessary 
for  large-scale  HP129Xe  production  suitable  for  applications  in 


Reserved  for  Publication  Footnotes 


[1]For  spin-1/2  particles  like  electrons  and  129Xe  nuclei,  polarization  is  defined  as:  P=(Na- 
Np)  /  (Na+Np),  where  Na  and  Np  refer  respectively  to  the  numbers  of  spins  in  the  a 
(mSi=+'\/2)  and  |3  (ms/  =-1/2)  magnetic  sublevels;  at  thermal  equilibrium,  PXe  is  or|ly 
~2.85*1Cr6  at  3  T  and  room  temperature. 
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Fig.  1.  (color  online)  The  XeNA  polarizer.  (A)Schematic  of  the  polarizer’s 
key  components  (the  liquid  N2  dewar  that  provides  gas  for  oven  heating 
and  cooling  is  not  shown).  The  optical  path  comprising  beam  expanding 
optics,  polarizing  beam-splitter  cube,  quarter-wave  plate,  and  heat  sinks  is 
represented  by  (“A/4”).  For  the  gas  cylinders,  “N”  and  “E”  designate  xenon  with 
naturally-abundant  129Xe  and  isotopically-enriched  129Xe,  respectively,  (B) 
Schematic  showing  the  principal  elements  of  the  optical  path. 


humans.  Here  we  present  results[2]  obtained  with  our  collab¬ 
oration's  first-generation  large-scale  (^1  L/hr)  129Xe  polarizer 
for  clinical,  pre-clinical,  and  materials  MRS/MRI  applications. 
Comprised  of  mostly  off-the-shelf  components,  our  automated, 
modular  polarizer  is  portable,  easy  to  use,  and  employs  an  ‘open- 
source’  design  that  is  readily  implementable  in  other  laboratories. 
Unlike  most  clinical-scale  Xe  polarizers,  our  hyperpolarizer  runs 
in  the  xenon-rich  gas  regime  (with  up  to  1800  Torr  in  a  500 
cc  optical  pumping  cell)  in  either  single-batch  or  stopped-flow 
mode,  in  part  negating  the  usual  requirement  to  cryo-collect 
the  HP  129 Xe — a  process  that  otherwise  increases  the  complexity 
of  the  device  and  can  lead  to  significant  losses  of  spin  polar¬ 
ization  during  HP129Xe  accumulation,  storage,  phase  transitions 
(33),  and  transfer.  Four  independent  methods  were  utilized  to 
measure  spin  polarization,  including:  in  situ  field-cycled  near-IR 
spectroscopy  (29)  (to  probe  Rb  electron  spin  polarization);  in 
situ  low-field  129Xe  NMR  (calibrated  with  thermal  XH  NMR);  gas 
transfer  to  47.5  mT  for  129Xe  NMR  (calibrated  with  thermal 13 C 
NMR);  and  gas  transfer  and  subsequent  129Xe  NMR/MRI  using 
a  clinical  3  T  scanner.  Excellent  agreement  was  observed  among 
these  different  approaches;  taken  together,  the  129Xe  polariza¬ 
tion  values  reported  here  represent,  by  a  significant  margin,  the 
highest  yet  achieved  at  such  high  Xe  densities — and  establish  the 
feasibility  of  attaining  near-unity  polarization  in  single  batches 
with  HP129Xe  quantities  sufficient  for  clinical  use.  Long  129Xe 
polarization  lifetimes  were  obtained  in  Tedlar  bags  that  are  well- 
explained  by  current  relaxation  models  (25),  and  that  are  more 


[2]Portions  of  this  work  were  presented  previously  in  Refs.  (31,  32). 
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Fig.  2.  (color  online)  In  situ  monitoring  of  global  Rb  electron  spin  polariza¬ 
tion.  Main  figure:  Examples  of  near-IR  spectra  of  the  pump  laser  transmitted 
through  the  cell  acquired  under  the  following  conditions:  room  temperature 
(“cold  cell”,  black  solid  lines);  57  °C  and  60=5.26  mT  (blue  solid  lines);  57  °C 
and  Bq=0  mT  (green  solid  Iinesy65  °C  and  fi0=5. 26  mT  (red  dashed  lines); 
and  65  °C  and  B0=0  mT  (orange  dashed  lines).  Three  scans  for  each  condition 
are  shown;  error  bars  for  (PRb>  values  were  determined  by  propagating  the 
standard  deviations  of  spectral  integral  values.  Insets: Cartoons  showing  the 
effects  of  depletion  pumping  of  electronic  spin  states  of  the  alkali  metal 
vapor  on  the  transmitted  laser  intensity  (assuming  OP  with  with  o+  circularly 
polarized  light;  nuclear  spin  levels  are  not  shown  for  simplicity).  When  the 
external  magnetic  field  B0  is  ‘off  (left),  optical  pumping  is  highly  inefficient, 
resulting  in  near-equal  populations  of  the  ground  electronic  spin  states 
(/77j=+1/2  &  -1/2)  and  a  relatively  high  density  of  absorbers  in  the  gas  phase. 
The  presence  of  an  applied  magnetic  field  provides  a  quantization  axis  and 
results  in  more  efficient  optical  pumping.  This  results  in  unequal  ground  state 
populations,  corresponding  to  high  electronic  spin  polarization;  the  reduced 
number  of  atoms  in  the  mj=-1/2  ground  state  gives  rise  to  increased  laser 
transmission  (right). 


than  sufficient  for  gas-phase  storage  or  accumulation  necessary 
for  some  applications. 

Methods 

XeNA  Polarizer.  Complete  details  of  the  design,  components,  assembly, 
and  operation  of  our  consortium’s  first-generation  ‘open-source’  automated 
129Xe  polarizer  (dubbed  “XeNA”  for  XEnon  polarizatioN  Automated)  will  be 
provided  elsewhere  in  a  full  report  (Nikolaou  eta!.,  manuscript  in  prepara¬ 
tion);  however,  brief  descriptions  of  key  aspects  are  provided  below.  XeNA 
is  self-contained  in  a  mobile  chassis  and  is  comprised  of  a  laser  system,  an 
electromagnetic  coil  assembly,  a  vacuum  /  gas-handling  manifold,  on-board 
optical  and  NMR  spectrometers  (and  other  sensors)  for  quality  assurance  (QA) 
and  feedback,  and  a  computer/automation  system  (Fig.  1A). 

The  optical  path  (Fig.  IB)  begins  with  a  200  W  laser  diode  array  (LDA, 
QPC  Lasers)  that  produces  frequency-narrowed  (30)  laser  output  at  the  Rb 
D-i  line  (~794.8  nm;  nominal  FWHM~0.27  nm)  The  laser  is  tunable  over  >1 
nm  by  adjusting  the  water-chiller  temperature  and/or  the  driving  current, 
allowing  the  laser  to  be  maintained  on  resonance  at  various  output  powers, 
as  well  as  off-resonant  SEOP  (29,  30,  34).  The  laser  beam  passes  through 
a  short  optical  fiber  preserving  linear  polarization  and  then  expands  to  2” 
diameter  before  collimation  and  entrance  into  a  2”  polarizing  beam-splitter 
(PBS)  cube.  The  s-polarized  beam  component  (<10%  of  incident  power) 
is  reflected  90°  and  discarded  into  a  beam  block;  the  main  (p- polarized) 
beam  is  transmitted  into  a  quarter-wave  (A/4)  plate,  which  renders  the  beam 
circularly-polarized.  This  beam  is  then  directed  into  the  oven,  which  contains 
a  cylindrical  glass  optical  pumping  cell  (2"  o.d.,  9.75"  long,  500  cc).  Prior 
to  installation  the  cell  is  coated  with  a  silanzing  agent  (SurfaSil,  Pierce)  to 
slow  129Xe  Ti  relaxation,  and  then  loaded  with  <200  mg  Rb  in  an  inert 
atmosphere.  SEOP  may  be  performed  once  the  cell  has  been  loaded  with  a 
variable  Xe/N2  gas  mixture.  The  N2  gas  helps  suppress  unwanted  re-emission 
of  non-polarized  light  from  the  electronically  excited  Rb  (1 5)  as  well  as  other 
undesirable  energetic  processes  (35);  the  N2  also  provides  some  additional 
collision-broadening  of  the  Rb  absorption  line  (36).  The  Teflon  oven  houses 
the  cell  and  provides  temperature  control  via  a  heated/cooled  gas  supply 
line  from  an  external  self-pressurized  liquid  N2  dewar.  The  laser  beam  enters 
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Fig.  3.  (color  online)  Determination  of  129Xe  polarization  at  low  and 
intermediate  magnetic  fields.  (A)HP  129Xe  NMR  spectrum  from  the  SEOP 
cell  containing  1545  Torr  Xe  (and  455  Torr  N2)  acquired  at  5.26  mT  (62.0 
kHz  129Xe  frequency;  1  scan  using  a  59  ps  RF  pulse  corresponding  to  a  30° 
tipping  angle — i.e.,  t30°=59  ps).  (fi)  Reference  NMR  spectrum  from  water  1 H 
spins.  111  M  1 H  concentration,  doped  with  5  mM  copper  sulfate,  thermally 
polarized  at  1 .46  mT  (62.0  kHz  1 H  frequency;  170,000  scans;  t30°=1  6  ps;  TR=0.3 
s;  Ph=5.0x10'9).  Inset:  1 H  FID  for  the  spectrum  in  the  main  figure;  SNR  was 
used  to  calculate  a  relative  error  bar  of  5.7%  for  PXe  values  calibrated  with 
this  signal.  (C)  HP  129Xe  NMR  spectrum  recorded  at  47.5  mT,  559  kHz  129Xe 
frequency,  1  scan,  using  t30°=18  psfrom  a  52  ml_  spherical  phantom  following 
transfer  of  some  of  the  cell  contents  (originally  300  Torr  Xe,  1700  Torr  N2)  to 
the  phantom.  (D)  Reference  13C  NMR  spectrum  from  170  mmol  of  sodium 
1-13C  acetate,  14  g  in  50  ml.  D20,  also  measured  at  47.5  mT  (508  kHz  13C 
frequency;  256  scans  using  t90°=54  ps;  TR=200  s;  Pc=  4.1  *10'8). 


time  (s) 

Fig.  4.  (color  online)  High-field  129Xe  measurements.  T ^  decay  of  normalized 
HP129Xe  NMR  signals  from  Xe  gas  transferred  to  Tedlar  bags  and  measured 
at  3  T  throughout  the  decay  period  (blue  triangles)  or  stored  at  ~1.5  mT 
(red  circles)  or  ~0  mT  (white  squares)  and  rapidly  transferred  to  and  from  3 
T  for  brief  acquisition  periods.  Loss  from  RF  excitation  pulsing  was  negligible, 
given  a  tipping  angle  of  <1°.  Inset.  Left:  129Xe  NMR  spectrum  from  HP129Xe  in 
a  Tedlar  bag  containing  ~800  cc  gas,  38%  Xe  by  volume  following  transport 
to  a  3  T  clinical  MRI  (34.09  MHz  129Xe  frequency;  one  scan  using  1 .4°  tipping- 
angle  RF  pulse).  Gas  compositions  for  this  spectrum  and  the  data  in  the  main 
figure  were  nearly  identical.  Right:  reference  NMR  spectrum  from  thermally 
polarized  129Xe  spins  in  the  3  L  spherical  phantom  containing  a  Xe/02  mixture 
(32  scans  with  90°  RF  pulses). 


and  exits  through  3”  anti-reflection  (AR)-coated  optical  flats.  The  beam  is 
retro-reflected  back  through  the  cell  via  a  3”  mirror  and  directed  into  the 
optics  assembly,  where  the  PBS  reflects  the  beam  into  a  second  beam  block. 
The  oven  and  cell  reside  in  a  homogeneous  magnetic  field  provided  by  a 
four-coil  electromagnet  assembly  (Acutran,  23.6"  ID;  nominal  fields  of  5.26 
or  1 .46  mT,  corresponding  to  62.0  kHz  resonance  frequency  for  129Xe  or 1 H, 
respectively).  The  entire  optics  train  and  magnet  assembly  is  contained  in  a 
light-tight  Alumalite  enclosure  (MiniTec). 

The  manifold  directs  gas  flow  from  high-pressure  cylinders  contain¬ 
ing  natural  abundance  xenon  (26.44%  129Xe),  isotopically  enriched  xenon 
(~86%  129Xe),  and  N2  gas,  providing  desired  mixtures  up  to  2000  Torr  total 
pressure.  All  experiments  used  naturally  abundant  129Xe  unless  otherwise 
stated.  Vacuum  is  provided  by  a  mechanical  rotary  vane  (‘rough’)  pump  and 
turbopump  station  (Edwards).  After  passing  through  02  getters,  the  gases 
are  loaded  into  the  OP  cell,  with  feedback  provided  by  digital  pressure 
gauges.  Following  SEOP,  the  HP129Xe  mixture  is  flowed  through  PFA  tubing 
and  a  Teflon  filter  (Millipore  Wafergard)  to  getter  any  residual  Rb  prior  to 
collection  in  a  Tedlar  bag.  The  bag  can  be  then  detached  and  transported 
elsewhere  for  experiments.  The  glass  valve  on  the  OP  cell  has  a  mechanically 
actuated  Teflon  stopcock.  The  other  manifold  valves  are  automated  pneu¬ 
matic  or  solenoid  valves,  where  Teflon  surfaces  are  used  for  all  valves  down¬ 
stream  of  the  OP  cell. 

In  situ  QA  for  the  SEOP  process  is  provided  by  near-IR  and  NMR  methods 
using  an  Ocean  Optics  HR4000  high-resolution  near-IR  optical  spectrometer 
and  a  Magritek  Kea2  low-field  NMR  spectrometer.  The  near-IR  spectrometer 
detects  via  a  fiber  probe  positioned  immediately  behind  the  3”  mirror.  The 
NMR  spectrometer  uses  a  home-built  NMR  probe/surface  coil.  Both  spectrom¬ 
eters  are  connected  to  the  laptop  computer  mounted  on  the  console.  User 
operation  of  XeNA  employs  a  GUI  designed  and  programmed  in  open-source 
software  (Processing.org),  which  drives  a  microcontroller  box  that  houses 
the  “brain”  of  the  polarizer  (an  Arduino  Mega  2560  REV3  microcontroller 
board),  solid  state  relays,  and  a  Burkert  manifold  enabling  pneumatic  valve 
operation.  The  micro-controller  is  responsible  for  control  of  all  solenoid 
and  pneumatic  Teflon  valves  as  well  as  communications  with  the  oven 
temperature  controller,  pressure  gauges,  laser  temperature  monitor,  and 
safety  interlocks.  Material  costs  for  the  hyperpolarizer  were  under  $125,000. 


Other  experimental  aspects.  Xenon  was  purchased  from  Nova  Gases.  MR 
spectra  and  images  were  obtained  at  47.5  mT  (37,  38)  or  3  T  field  using  a 
Magritek  Kea2  with  permanent  magnet  (559  kHz  129Xe  frequency  and  508 
kHz  13C  frequency)  or  a  Siemens  3  T  clinical  MRI  scanner  (34.09  MHz  129Xe 
frequency),  respectively.  Elemental  analysis  of  Tedlar  bag  contents  (following 
SEOP  runs)  was  performed  by  Element  One  (Wilmington,  NC)  to  ensure  that 
<5.0  ng  Rb  was  present  in  the  bag  after  transfer. 

Results 

An  estimate  of  the  spatial  average  of  the  Rb  electron  spin 
polarization,  (Pro),  can  be  obtained  from  measurement  of  the 
transmitted  light  of  the  pump  laser  while  the  magnetic  field  is 
cycled  (29)  (Fig.  2).  Spectra  are  recorded  from  the  laser  beam 
transmitted  through  the  cell  in  order  to  calculate  absorbance 
values: 


where  A  and  T  respectively  denote  the  absorbance  and  trans¬ 
mittance,  and  /hot  and  /coid  are  the  integrated  intensities  of  the 
transmitted  laser  spectra  obtained  when  a  gas-loaded  cell  is 
respectively  ‘hot’  and  ‘cold’  (i.e.,  with  and  without  significant 
Rb  vapor  present).  With  some  simple  assumptions,  absorbance 
values  may  then  be  compared  under  the  same  conditions — except 
with  or  without  the  SEOP  magnetic  field  (Bo)  on — in  order  to 
calculate  an  estimate  of  (Pro),  according  to  a  simple  relation 
derived  from  Beer’s  Law  (29): 
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Fig.  5.  (color  online)  129Xe  nuclear  spin  polarization  values  measured  at 
5.26  mT,  47.5  mT,  and/or  3  T,  plotted  versus  xenon  partial  pressure  prior 
to  SEOP.  Labels  ‘b.t.’  and  ‘a.t.’  respectively  refer  to  measurements  obtained 
from  Xe  gas  remaining  within  the  SEOP  cell  before  and  after  some  of 
the  gas  was  transferred  to  another  container.  Error  bars  are  determined 
from  the  uncertainties  in  the  spectral  integral  values  obtained  from  the 
respective  thermally-polarized  reference  samples.  The  value  at  725  Torr  was 
obtained  with  82%-enriched  129Xe.  Inset.  2D  FLASH  129Xe  projection  image 
of  hyperpolarized  xenon  gas  in  a  Tedlar  bag  following  transport  to  the  3  T 
MRI  (FOV:  40  cm;  slice  thickness:  1 .7  cm;  matrix  size:  80x80;  TE/TR:  2.6/5. 1  ms; 
RF  tipping  angle  ~ 4°;  and  SNR  ~40:1). 


whereto  is  the  absorbance  when^o =0.  Examples  of  transmission 
spectra  are  shown  in  Fig.  2  for  a  cell  containing  495  and  1300 
Torr  of  Xe  and  N2,  respectively.  Multiple  near-IR  spectra  were 
obtained  in  rapid  succession  under  a  given  set  of  conditions,  and 
then  averaged  and  integrated.  One  set  of  scans  was  obtained  at 
57  °C,  roughly  10  minutes  into  the  SEOP  process;  the  corre¬ 
sponding  values  of  (Pm,)  and  Pxe  measured  with  low-field  NMR 
(see  below)  were  43±9%  and  48.3±2.8%,  respectively. [3]  As  part 
of  the  SEOP  optimization  process,  the  temperature  of  the  cell 
was  raised  to  65  °C  and  after  several  minutes  the  sequence  was 
repeated.  Improved  values  of  (Pm>)  andPxe  of  71±6%  and57±3%, 
respectively,  were  obtained. 

The  Xe  polarization  process  was  monitored  directly  via  NMR 
under  three  different  sets  of  conditions  (Figs.  3&4).  First,  low- 
field  129Xe  NMR  can  be  used  to  measure  Pxe  within  the  OP 
cell — during  or  after  the  SEOP  process.  Fig.  3A  shows  an  example 
of  a  HP129Xe  NMR  spectrum  obtained  at  5.25  mT  from  a  cell 
containing  1545  Torr  Xe  following  SEOP  and  cool-down.  Low- 
field  NMR  is  calibrated  using  a  thermal  1H  reference  signal  ob¬ 
tained  separately  using  the  same  NMR  circuit  at  the  same  Larmor 
frequency  of  62  kHz,  and  tuning/matching  conditions  detecting 
XH  at  1.46  mT  from  an  identical  cell  filled  with  5  mM  Q1SO4 
doped  water  (Fig.  3B).  As  part  of  our  initial  tests  of  the  efficiency 
of  the  Xe-transfer  process,  Pxe  measurement  was  also  performed 
at  an  intermediate  field  of  47.5  mT  located  approximately  2  m 
away  from  the  polarizer.  Fig.  3C  shows  an  example  of  a  HP  129Xe 
NMR  spectrum  obtained  at  this  field  from  a  portion  of  HP129Xe 
within  a  cell  containing  300  Torr  Xe  following  automated  transfer 


to  a  ~50  cc  polypropylene  spherical  phantom  connected  to  the 
polarizer  via  %”  OD  PFA  tubing  and  located  within  the  field 
of  the  permanent  magnet.  129Xe  polarization  (here,  68.7±3.4%) 
was  calibrated  using  a  thermal  13  C  reference  signal  obtained 
from  170  mmol  of  sodium  1-13C  acetate  in  99.8%  D20  located 
in  an  identical  phantom,  Fig.  3D.  Comparison  with  in  situ  5.26 
mT  129Xe  NMR  polarimetry  taken  from  the  OP  cell  during  the 
same  experiment  before  transferring  the  gas  to  47.5  mT  magnet 
(68.5±3.9%)  indicated  no  polarization  loss  within  error,  and  thus 
a  highly  efficient  HP129Xe  transport  process. 

The  hyperpolarizer  was  delivered  to  Brigham  &  Women’s 
Hospital  in  Boston,  MA  and  installed  adjacent  to  a  clinical  MRI 
suite  in  February  2012.  Optimization  of  the  cell  cool-down  and 
automated  Xe  transfer  process  allowed  Xe  transfer  to  Tedlar 
bags  via  expansion  and  subsequent  bag  transport  to  a  3  T  clinical 
MRI.  The  inset  of  Fig.  4  shows  an  example  of  a  high-field  HP 
129Xe  NMR  spectrum  from  a  Tedlar  bag  containing  ~800  cc  gas 
(38%  Xe  by  volume,  given  an  initial  cell  loading  of  760  Torr  Xe  / 
1240  Torr  N2).  The  Pxe  value  is  calibrated  using  a  thermal  129Xe 
NMR  signal  from  a  3  L  spherical  reference  sample  containing  4 
atm  Xe  with  86%  129Xe  enrichment  and  2  atm  02  with  a  T\  of 
0.99  s.  Again,  little  loss  of  polarization  was  suffered  throughout 
the  transfer  /  transport  process,  given  values  of  42.8±2.4%  and 
41±2%  for  Pxe  measured  in  the  cell  prior  to  transfer  and  at  3 
T  in  the  Tedlar  bag,  respectively.  The  latter  %Pxe  corresponds 
to  a  ~  144,000-fold  polarization  enhancement  over  the  thermal 
equilibrium  value.*  The  strong  NMR  signals  obtained  from  the 
HP  129Xe  gas  allowed  the  polarization  decay  to  be  monitored 
under  a  variety  of  conditions.  Figure  4  shows  129Xe  T 1  relaxation 
data  for  HP129Xe  in  Tedlar  bags  recorded  at  3  T,  but  stored  at 
different  field  strengths:  following  exponential  fits,  T 1  decay  times 
of  38±12  min  and  5.9±0.4  hr  were  observed  for  129Xe  gas  in  Tedlar 
bags  at  1.5  mT  and  3  T,  respectively.  Such  long  T\  values  bode 
well  for  a  variety  of  applications  where  gas-phase  storage  and/or 
accumulation  of  HP129Xe  is  needed  (25). 

129Xe  polarization  values  for  various  Xe  densities,  operating 
conditions,  and  measurement  magnetic  fields  are  summarized  in 
Fig.  5.  The  data  exhibit  good  run-to-run  and  method-to-method 
agreement  across  the  various  measurements  for  the  given  Xe  cell 
loadings  (^300-1600  Torr),  despite  the  fact  that  the  Px&  values 
were  measured  at  different  fields,  with  different  NMR  acqui¬ 
sition  methodologies,  and  different  thermal  reference  samples. 
The  data  include  values  for  in-cell  Pxe  values  of  90.9±5.2%, 
57.1±3.3%,  50.1±2.9%,  and  33.4±1.9%  measured  for  Xe  load¬ 
ings  of  300,  495,  765,  and  1570  Torr,  respectively;  Px&  values  of 
41±1.6%  and  28±1.1%  with  ~760  and  ~1545  Torr  Xe  loadings 
were  obtained  following  transfer  to  Tedlar  bags  and  subsequent 
transport  to  the  3  T  scanner,  permitting  imaging  with  high  SNR 
(~40)  to  be  demonstrated,  Fig.  5  inset.  Moreover,  these  post¬ 
transfer  values  are  ~2.7  times  greater  than  previous  bests  in  Pxe 
achieved  near  such  high  Xe  densities;  when  combined  with  the 
increase  in  cell  volume,  these  results  constitute  a  ~  18-fold  im¬ 
provement  in  HP129Xe  production  over  the  previous  best  results 
obtained  at  such  high  in-cell  Xe  densities  (23). 

Discussion 

The  overall  performance  of  the  polarizer  is  the  result  of  a  num¬ 
ber  of  aspects  of  the  SEOP  process  that  are  both  fundamental 
and  technical  in  nature  (30).  SEOP  can  be  treated  as  a  simple 
relaxation  process  (24),  which  at  steady  state,  simplifies  to  give 
(16): 


[3]PXe  cannot  be  higher  than  (PRb)  (see  Eq.  (3));  however,  the  difference  between  these 

values  is  less  than  the  uncertainty.  •(/>») 

/s£  1  ,  (3) 
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where  Yse  is  the  Rb/Xe  spin  exchange  rate  and  r xe  is  the  129Xe 
nuclear  spin  destruction  rate  (= l/7"i);  thus,  it  is  convenient  to 
categorize  aspects  of  the  hyperpolarizer  and  its  operation  as 
factors  that  respectively  help  optimize  ( Pro ),  Yse,  and/or  Txe- 
First,  the  most  important  factor  in  maximizing  global  (PRb)  is  the 
cell  illumination  by  the  laser.  At  a  given  position  (r)  within  the 
cell,  the  local  PR b  is  given  by  (16): 


Yop(r) 

r0p(r)+YsD 


.(4) 


Where  Yop(^)  is  the  local  optical  pumping  rate  (given  by  the 
integrated  product  of  the  laser  flux  and  the  Rb  absorption  cross- 
section  (36)),  and  Tsd  is  the  Rb  electronic  spin  destruction 
rate,  which  is  dominated  by  non-spin-conserving  collisions  with 
Xe  (specifically,  via  the  spin-rotation  interaction  (26,  28))  and 
is  high  under  our  conditions  because  of  its  proportionality  to 
[Xe].  Correspondingly,  the  laser  power  employed  here  is  also 
high — about  170  W  for  most  experiments,  mostly  resonant  or 
near-resonant  with  the  Rb  Di  line.  In  fact,  the  resonant  photon 
flux  is  roughly  six-fold  greater  than  in  Ref.  (23),  but  this  increase 
explains  only  part  of  the  improvement  in  HP129Xe  performance. 
Care  was  also  taken  with  the  optical  path  to  ensure  complete 
illumination  throughout  the  optical  pumping  cell  with  near-unity 
laser  polarization,  with  no  dark  regions  (either  near  the  edges 
or  from  longitudinal  optical  nodes),  and  with  additional  photon 
flux  provided  by  retro-reflection.  Control  of  the  cell  temperature 
is  also  key  for  SEOP  optimization.  Previously  it  was  shown  that 
there  can  be  an  inverse  relationship  between  Xe  density  and 
the  optimal  temperature  for  SEOP,  where  higher  [Xe]  favors 
lower  cell  temperatures  (23, 30) — an  effect  that  may  be  explained 
in  part  by  the  need  to  maintain  uniform  illumination  through¬ 
out  the  cell.  Although  Yse  is  proportional  to  the  Rb  density 
(39) — ostensibly  favoring  higher  cell  temperatures — maintaining 
a  high  “photon-to-Rb”  ratio  (in  part  by  limiting  [Rb])  can  be 
particularly  important  for  preserving  global  (Pro)  as  Xe  density 
is  increased,  and  Xe-induced  Rb  spin-destruction  becomes  domi¬ 
nant.^]  For  example,  ongoing  simulations  (40)  predict  a  high  F sd 
value  of  ~  134,300  s'1  for  the  experimental  conditions  of  the  65 
°C  data  in  Fig.  2,  dominated  by  Xe  collisions  (26-28).  However, 
the  high  resonant  photon  flux  at  the  front  of  the  cell  should  give 
nearly  a  9-fold  greater  optical  pumping  rate  (yop~1.16x105 6  s'1); 
using  Eq.  (4),  such  values  for  F sd  and  Yop  would  correspond  to 
Pro ^89.6%.  Maintaining  a  relatively  low  Rb  density  allows  the 
laser  to  penetrate  efficiently  to  ensure  high  photon  flux  and  Rb 
polarization  throughout  the  cell,  consistent  with  the  experimental 
value  ((P Rb)  =71±6%).  Thus,  particularly  in  the  limit  where 
r Xe  is  negligible,  maximizing  (PRb)  is  more  important  than  Yse, 
once  Yse  is  sufficiently  high.  Indeed,  our  high  Pxe  values  were 
enabled  in  part  by  our  ability  to  achieve  long  in-cell  129Xe  T\ 
(=l/rxe)  decay  times.  As  one  example,  following  SEOP  with  a 
cell  containing  495  Torr  Xe  and  cool-down  to  33  °C  (where  Rb 
should  be  condensed),  a  cell  T\  of  1.9±0.6  hr  was  measured  at 
5.26  mT.  Given  a  typical  build-up  time  constant  measured  for 
one  experiment  of  (YsE  +  Fxe)"1  =  8.5  min,  this  in-cell  129Xe  T\ 
translates  to  a  spin-exchange  rate  of  Yse^T.8x10'3  s'1 — roughly  an 
order  of  magnitude  greater  than  Txe.  Thus,  according  to  Eq.  (3) 
Pxe  should  closely  approach  (PRb)  under  these  conditions,  again 
in  good  agreement  with  the  Fig.  2  experiments. 

In  addition  to  long  HP129Xe  relaxation  times  measured  in  OP 
cells,  very  slow  polarization  decay  was  also  observed  following 


[4]An  additional  contribution  to  the  interplay  of  Xe  density  and  SEOP  temperature  may 
be  poor  energy  dissipation  caused  by  reduced  thermal  conductivity  of  the  gas  for  Xe- 
rich  mixtures — a  possibility  we  are  currently  studying  with  in  situ  Raman  spectroscopy 
(Whiting,  Newton  etal.,  manuscripts  in  preparation). 


transfer  to  Tedlar  bags.  Generally,  129Xe  relaxation  rates  are 
determined  by  several  potential  contributions  (adapted  from  Ref. 
(25)): 


r*  =  ri+rs+r02+K  ,(5) 

where  T*  is  the  ‘intrinsic’ — and  essentially 

unavoidable — contribution  from  spin-rotation  interactions 
of  transient  and  persistent  Xe-Xe  dimers  (Vi  =  Vt  +  Vp  (41)),  Vg 
results  from  diffusion  through  field  gradients,  T 02  results  from 
dipolar  interactions  with  residual  paramagnetic  O2  in  the  gas 
mixture,  and  (Tw  =  Vw>e-  +  r^Ar)  includes  contributions  from 
Xe/wall  collisions  involving  interactions  with  surface  spins  of 
unpaired  electrons  and  nuclei,  respectively.  Vw>e.  can  be  neglected 
for  Tedlar  (polyvinyl  fluoride);  moreover  Tg  can  be  neglected  for 
the  3  T  measurements  given  the  homogeneous  field.  Saam  and 
co-workers  recently  provided  a  semi-empirical  formula  for  T; 
that  when  adapted  for  our  conditions,  is  given  by  (25): 


r.  = 


[Xe] 


56.1  hr  4.59  hr 


[l  +  (3.65xlO-;)Si]l  +  0.51^j 


(6) 


where  the  two  additive  terms  respectively  correspond  to  T t  and 

IT^,  Bo  il  in  Tesla,  the  gas  densities  are  in  amagat[5],  and  0.51  is 
a  factor  that  takes  into  account  the  differential  break-up  rate  of 
persistent  Xe  dimers  from  collisions  with  N2.  Using  values  for 
[Xe]  and  [N2]  of  0.35  and  0.58  amg,  Eq.  (6)  predicts  a  limiting 
value  of  l/\~xe  =7.8  hr  (with  10%  uncertainty),  in  relatively  good 
agreement  with  our  experimental  value  of  5.9±0.4  hr  at  3  T, 
and  indicating  effective  suppression  of  other  relaxation  pathways. 
Given  the  expected  absence  of  paramagnetic  wall  sites,  if  we 
assume  that  (a)  relaxation  from  wall  collisions  involving  surface 
nuclear  spins  (42)  (e.g.  and  19F)  is  essentially  quenched  at 
high  field  and  (b)  that  the  remaining  contribution  to  Vxe  is  from 
residual  O2  (with  a  relaxivity  of  ^0.4  Hz/amg  (43)),  then  the 
differential  relaxation  rate  would  correspond  to  an  upper  limit  for 
the  O2  partial  pressure  of  ~2X10'2  Torr — a  reasonable  amount. 
While  Vp  is  the  limiting  factor  for  our  high-field  129 Xe  relaxation, 
the  accelerated  decay  at  1.5  mT  (Ti=38±12  min)  is  consistent 
with  VW)x  providing  the  dominant  mechanism.  Driehuys  et  al.  (42) 
showed  that  at  low  field  (~0-10  mT),  dipolar  relaxation  with 
XH  ‘wall’  spins  was  the  primary  contributor  to  129Xe  polarization 
decay  in  cells  with  organosilane  coatings,  with  Ti<50  min  at 
1.5  mT  and  297  K.  Although  Tedlar  may  lack  the  nanoscale  Xe 
permeability  of  SurfaSil  coatings  (42),  it  does  possess  a  high  sur¬ 
face  density  of  and  19F  spins  and  likely  provides  an  attractive 
surface  for  transient  Xe  adsorption.  More  detailed  studies  of 
129Xe  relaxation  in  Tedlar  containers  will  be  the  subject  of  future 
efforts. 

On  a  technical  level,  SEOP  efficiency  is  aided  by  the  presence 
of  in  situ  monitoring  of  the  Rb  spectral  absorption  and  129Xe 
NMR  signals,  which  provide  immediate  feedback  to  allow  opti¬ 
mization  of  the  SEOP  process  in  real  time.  Afterwards,  additional 
129Xe  polarization  losses  are  mitigated  by  (a)  rapid  cool-down  of 
the  cell  while  reducing  the  resonant  laser  power  in  kind  (to  further 
decrease  the  heat  load,  while  still  providing  sufficient  illumination 
to  avoid  “dark”  Rb — a  procedure  that  allows  only  a  few  percent 
loss  of  Pxe),  as  well  as  ( b )  suppressing  exposure  to  O2  and  other 
paramagnetic  materials  throughout  the  gas  lines,  (c)  moving  the 
HP129Xe  quickly  into  the  sample  or  transport  vessel,  and  (< d ) 
avoiding  Xe  phase  transitions  and  other  relaxation-susceptible 
portions  of  the  Xe  phase  diagram  (33).  Indeed,  Xe  cryo-collection 


[5]  1  amagat  (amg)  =  2.6873  *  1019  cm3,  equal  to  the  density  of  an  ideal  gas  at  0  °C  and 

760  Torr. 
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/  sublimation  prior  to  transfer  to  the  sample  is  optional  for  XeNA, 
but  it  is  normally  a  requirement  with  other  polarizer  designs 
operating  with  standard  “xenon-lean”  gas  mixtures.  We  should 
add  that  not  cryo-collecting  the  Xe  does  result  in  dilution  of  the 
HP129Xe  with  N2  gas  and  leaves  behind  significant  fractions  of 
the  SEOP  mixture  in  the  cell  following  the  expansion  process; 
however,  both  of  these  issues  can  be  mitigated  respectively  by 
using  Xe-rich  mixtures  and  by  the  addition  of  a  large  automated 
gas  piston  (16)  where  the  cell  contents  can  be  expanded  into  a 
much  larger  volume  prior  to  transfer  to  the  sample  or  transport 
vessel. 

Conclusion 

We  have  presented  results  from  a  clinical-scale  (~1  L/hr) 
hyperpolarizer  that  produces  batches  of  spin-polarized  129Xe 
sufficient  for  clinical  use.  The  polarizer’s  ‘open-source’  design 
and  automated  operation  should  facilitate  implementation  of 
HP129Xe  technology  into  other  labs,  including  clinical  settings. 
Good  agreement  was  obtained  among  the  four  independent 
methods  for  characterizing  the  spin  polarization.  The  ability  to 
achieve  high  Px&  values  at  high  Xe  densities,  combined  with 
stopped-flow  operation,  negates  the  usual  requirement  of  Xe 
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cryo-accumulation  and  storage,  and  opens  a  door  to  greatly  im¬ 
proved  polarization  efficiency  for  quadrupolar  isotopes  (e.g. 83 Kr 
and  131Xe  (44,  45)).  Further  improvements  in  design  and  gas¬ 
handling  efficiency  will  be  manifested  in  our  2nd -generation  hy¬ 
perpolarizer  (now  under  development);  other  possible  improve¬ 
ments  in  efficiency  for  large  scale  hyperpolarization  of  the  heavy 
noble  gases  may  exploit  Cs  (e.g.  Ref.  (46))  or  Cs/Rb  hybrid  SEOP 
schemes.  These  results,  combined  with  ultra-long  polarization 
lifetimes  in  the  gas  phase,  bode  well  for  in  vivo  experiments. 
Indeed,  recent  FDA  IND  and  IRB  regulatory  approval  for  the 
polarizer  will  enable  a  host  of  planned  experiments  that  include 
a  variety  of  clinical  applications  with  human  subjects. 
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